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(N ' Abstract 
> 
O 

This is the first in a series of two papers. In this first part, we use the Schwinger-Keldysh formalism to 
{Sj ' derive semiclassical Boltzmann transport equations, accurate to order h, for massive chiral fermions, 
scalar particles, and for the corresponding CP-conjugate states. Our considerations include com- 
O ' plex mass terms and mixing fermion and scalar fields, such that CP-violation is naturally included, 
I rendering the equations particularly suitable for studies of baryogenesis at a first order electroweak 
phase transition. We provide a quantitative criterion for when the reduction to the diagonal kinetic 
^| equations in the mass eigenbasis is justified, leading to a quasiparticle picture even in the case of 
^ ; mixing scalar or fermionic particles. Within the approximations we make, it is possible to first study 
the Boltzmann equations without the collision term. In a second paper |lj we discuss the collision 
^ ■ terms and reduce the Boltzmann equations to fluid equations. 
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1. INTRODUCTION 



The question of how Boltzmann transport equations emerge from quantum field theory is an old 
problem 0,H,0]- It seems that a consensus has been reached, according to which a suitable (loop) 
truncation of the self-energies in the Kadanoff-Baym equations, supplied by an on-shell reduction, 
leads to Boltzmann transport equations for the dynamics of quasiparticles in phase space. While 
procedures leading to classical Boltzmann equations - which formally correspond to taking the limit 
h — > of the Kadanoff-Baym equations - seem to be quite clear and well established [1, , no system- 
atic treatment has been pursued beyond the classical approximation, which would lead to a kinetic 
description valid to linear order in an h expansion (a notable exception is a study of interference 
between the tree- level and one-loop decay rates of relevance for grand unified baryogenesis [7].) The 
inclusion of h corrections is absolutely necessary for any treatment of baryogenesis, which represents 
the primary impetus for the work presented here, since CP-violating effects vanish in the limit when 
h — > 0. More concretely, one deals with particle transport in high temperature relativistic plasmas 
with Minimal Standard Model particle content, often extended to include exotic, as-of-yet undiscov- 
ered new particles, e.g. various supersymmetric partners. While h expansions are often identified 
with loop expansions in field theory, here we associate h primarily to a gradient expansion (which 
is a field-theoretical generalization of the WKB expansion) with respect to space-time variations of 
background fields. The relevant examples of background fields for baryogenesis are the mass terms 
generated by the (Higgs) order parameter at a first order electroweak phase transition; other exam- 
ples include classical (possibly stochastic) background gauge fields. One important hurdle, which 
must be overcome in the development of any formalism that aspires to accurately describe transport 
in plasmas beyond the classical level, is the establishment of a quasiparticle picture. This entails a 
demonstration of the existence and the identification of the relevant quasiparticle excitations in the 
plasma, which can be used to encode the relevant dynamical information. This is the problem we 
solve here. Apart from baryogenesis, potential applications of the formalism developed in this work 
include propagation of neutrinos in dissipative media, transport of electrons and holes in quantum 
wires and quantum semiconductor devices, dynamics of Bose-Einstein condensates, etc. 

The standard approach to electroweak baryogenesis requires a study of the generation and trans- 
port of CP-violating flows arising from interactions of fermions with expanding phase transition 
fronts of a first order electroweak phase transition. The most prominent theoretical problem in such 
scenarios over the past few years has been to find a systematic derivation of the appropriate trans- 
port equations, including the CP-violating sources. In particular, there has been controversy in the 
literature regarding what are the dominant sources appearing in the transport equations and how to 
compute them. This work is an attempt to resolve this daunting question. 

It is known that a strong first order transition is a necessary requirement for electroweak baryogene- 
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While absent in the Minimal Standard Model 
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15, 



can be realized in its extensions, which include two Higgs d oublet models 



symmetric Standard ModelfMSSM) Q, 0, H Q HQ H H E 



extensions (NMSSM) 



It Jl7J, a strong transition 
J, the Minimal Super- 

|, and its nonminimal 

I32] ]. In the NMSSM, not only the possibility of electroweak baryoge- 
nesis is given in wide regions of the parameter space, but it can also have stronger CP-violation, 
increasing the amount of the actually produced baryon asymmetry [3^ . 
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Considering the derivation of the relevant transport equations, there are essentially two approaches 
in literature. The first approach, which we shall refer to as the semiclassical force mechanism 0,0, 
is based on the observation that, when treated in the WKB approximation, a CP- violating quadratic 
part of the fermionic Lagrangean exhibits a CP-violating shift at first order in gradients in the 
dispersion relation, which manifests itself as a CP-violating semiclassical force in the Boltzmann 
equation. This force then appears in the fluid transport equations, or equivalently diffusion equations 
for the relevant particle densities. In the original work |34| |35( baryogenesis from scatterin g of top 
quarks and tau leptons in two Higgs doublet models was considered. In the subsequent work [3i| the 
semiclassical force formalism was extended to include the case of fermion mixing, which is relevant 
for example for the charging induced baryogenesis in the Minimal Supersymmetric Standard Model 



(MSSM) (3J, |3jJ |3J, |39|, |4(| and in its nonminimal extensions (NMSSM) (41 
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The second approach, commonly referred to as spontaneous baryogenesis [43], is based on the obser- 
vation that in the presence of CP-violation the (fermion) hypercharge currents are not conserved. 
As a consequence, the energy levels for CP-conjugate fermionic states are mutually shifted. In the 
presence of scatterings and transport, this shift then leads to different populations for chiral fermions 
and the corresponding antifermions, thus sourcing baryogenesis. It was then pointed out in Refs. 0], 
and that the hypercharge current, which sources spontaneous baryogenesis, is suppressed 
by the square of the mass. (This must be so simply because the hypercharge current is conserved 
in the limit when the mass of particles vanishes, so that the spontaneous source must vanish, too.) 
Mostly in order to incorporate the possibility of fermionic mixing, the spontaneous baryogenesis 
mechanism was subsequently developed by several groups, which adopted the idea to baryogenesis 
studies mediated by the stops and charginos of the MSSM 00, 00,0,0 and NMSSM 0. 

From the current literature, it is in fact not clear whether the semiclassical force and spontaneous 
sources represent identical or different sources, or if there perhaps has been a conceptual problem in 
calculating (one or both of) the sources. The situation was made even more controversial by a recent 
work [51], in which it was claimed that the spontaneous baryogenesis and semiclassical force sources 
correspond formally to the same source in fluid transport equations, arguing that the method used 
for computing the semiclassical force source was incorrect. 
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This is to be compared with the methodology advocated in the work based on the WKB method 
and more recently in 0000 based on the Schwinger-Keldysh formalism |2| , or more precisely 
on a gradient expansion of the Kadanoff-Baym equations [3] in a weak coupling regime. In order to 
clarify the point of disagreement, we first extend the approach developed in |52j, |^3J. We consider 
the dynamics of non-equilibrium Wigner functions for massive chiral fermions which couple to a CP- 
violating Higgs condensate of advancing phase fronts to first order in a gradient expansion. Assuming 
that the relevant mass parameter depends only on one spatial coordinate, which models reasonably 
well the limit of large almost planar bubble walls of a strong first order phase transition, implies a 
conserved spin, and the constraint equation for the Wigner function is solved by a spectral on-shell 
ansatz, thus proving the validity of the quasiparticle picture of the plasma. This result holds both in 
the case of one fermion and in the case of mixing fermions. 

An important difference, regarding the source calculation in spontaneous baryogenesis and semi- 
classical force baryogenesis, is the choice of the basis in flavor space, which leads to qualitative 
differences in the resulting baryon production. In studies of chargino-mediated baryogenesis flavor 
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mixing is of crucial importance. While the proponents of spontaneous baryogenesis argue that the 
weak interaction (flavor) basis (of charginos) is the right basis to work in, since it is favored by 
the interactions [51], the semiclassical force camp contends that the mass eigenbasis has apparent 
advantages. Of course, any physical quan tity must be basis independent, but only when one includes 
flavor mixing in transport equations [56|, which - as of yet - has not been done. We strongly favor 
the mass eigenbasis: as we prove in section 15.21 when working in the mass eigenbasis and to order h 
accuracy, the diagonal densities decouple from the off-diagonals. The decoupling works uniquely in 
the mass eigenbasis. The proof is valid provided the mass eigenvalues are not nearly degenerate, that 
is when hk-d -C S(m^) is satisfied, where S(rrv^) signifies the split in the mass eigenvalues. While this 
strongly indicates that the mass eigenbasis is singled out, the definite resolution of the controversy 
awaits a basis independent treatment. 

The paper is organized as follows. 

In section 121 we review the Schwinger-Keldysh formalism, suitable for out-of-equilibrium dynamics 
of field theoretical correlators, and present a derivation of the relevant Kadanoff-Baym equations in 
Wigner space for our model Lagrangean of chiral fermions coupled to scalars. We include the case 
of flavor mixing in both the fermionic and the scalar sector. 

In section El we study an on-shell reduction of the Kadanoff-Baym equations for a single scalar field. 
We show that it is consistent to include the (real part of the) self-energy and the collision term to 
both the constraint and kinetic equation, provided one truncates the self-energy and the collision 
term at the same order in the coupling constant. This resolves one of the important - and up-to- 
now unanswered - questions of Boltzmann-like kinetics of relativistic systems. Even though we do 
not attempt to generalize our result to the kinetics of fermions, it is quite plausible that the same 
conclusion can be reached. In the subsequent sections we assume that this is the case. 

A proof that no source appears in the flow term of the kinetic equations for mixing scalars at order 
h for planar walls is presented in section 0] We complete the proof originally presented in 0|. Next, 
we discuss scalar kinetic equations for the CP-conjugate scalar densities, which are then used in 
Paper II [1] to identify the CP-violating source in the scalar collision term. 

Section is devoted to the tree-level kinetics of mixing fermions in the presence of a CP-violating, 
complex mass matrix. An important application includes chargino and neutralino mediated baryo- 
genesis in supersymmetric extensions of the Minimal Standard Model. We begin by identifying a 
conserved quantity, which for planar bubble walls and in the wall frame corresponds to the spin in 
the direction ft oc (k x k z , k y k z , k\ — k\ — ky) , which does not correspond to helicity, which has been 
erroneously identified as the conserved quantity in most of the literature on semiclassical force baryo- 
genesis. Then we construct Wigner functions which are block diagonal in the conserved spin. Just 
like in the scalar case, we prove that to order h the diagonal and off-diagonal densities decouple in the 
mass eigenbasis (we fill in a gap in the original derivation in Ref. provided the mass eigenvalues 
are not nearly degenerate, hk ■ d <C 5(\rrid\ 2 )- The final part of the section contains the derivation of 
Boltzmann-like transport equations for the fermion distribution functions, with particular emphasis 
on the kinetics of CP-violating densities. The relevant CP-violating sources are identified, some of 
which are unaccounted for in the existing literature. For example, we show that a CP-even deviation 
from equilibrium in the distribution function can contribute as a CP-violating source at order H. 

In Paper II [1] we address the collision terms of the Boltzmann equations, which contain further 
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CP-violating sources. We furthermore simplify the Boltzmann equations to a set of fluid equations, 
based on which one can study the dynamics of CP-violating currents. 
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CL 1 

FIG. 1: The complex time contour for the Schwinger-Keldysh non-equilibrium formalism. 

2. FROM THE 2PI EFFECTIVE ACTION TO KINETIC THEORY 
2.1. Green Functions 

The in- out- formalism of quantum field theory is well suited for the description of particle scattering 
experiments, in which the system is prepared to be in a definite in-state at t — ► — oo and the question 
is what is the probability amplitude for the system to be found in a definite out-state at t — > +oo. 
In statistical physics, however, it is often of interest to follow the temporal evolution of a system. 
Starting with definite initial conditions at some time t = to, we ask for the expectation values of 
physical quantities at finite times. A theoretical framework for such problems was first suggested by 
Schwinger in 1961 Q and then developed further by Keldysh j^J and others. An extension of field 
theory capable of dealing with non-equilibrium problems is obtained by defining the time arguments 
of all quantities on a path C that leads from the initial time t Q to t and then back to t , as illustrated in 
figure All integrals and derivatives have then to be performed along that path, and the usual time 
ordering becomes time ordering Tq along C. This formalism is called Schwinger-Keldysh formalism 
or Closed Time Path (CTP) formalism p|. The definitions of the scalar and the fermionic Green 
functions are 

iA(u,v) = (n\T C (j)(u)^(v)\n} (2.1) 
iS af) (u,v) = (n\T c iP a (u)Mv)\n) , (2.2) 

where \ denotes the physical state of the system, and and if) are the scalar and fermionic fields, 
respectively. The contour C can now be split into a C + branch from to to t and a C_ branch from t 
to to, as shown in figure ^ An important simplification, which will help us to perform the Wigner 
transform, instrumental for gradient expansion, is an extension of the integration path to to ~~ > — oo 
and t — > oo. This simplification is suitable for problems of our interest, which comprise plasmas close 
to chemical and thermal equilibrium with efficient equilibration processes, as it is indeed the case with 
the electroweak plasma. In this case the influence of initial conditions can be neglected. Denoting 
the branch on which the time argument lies by an index a = ±, we can rewrite the formalism using 
ordinary time arguments. We then have 



/ d A u — > > a / d 4 u 

Jc V loo 

8c{u — v) -> a5 ab 5(u-v) 

iA(u,v) -> iA ab (u,v) 



iS(u,v) -> iS ab (u,v). (2.3) 

The additional factors of a in the integral and in the 5-function appear since the C_ branch runs 
backward in time. 
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In this Keldysh formulation one usually defines the following four scalar Green functions: 

iA ++ (u,v) = zA*(u,v) = (n\T[<t>{u)<t>\v)]\n) 
iA + -(u,v) = iA < (u,v) = (n\^(v)(f)(u)\n) 
iA- + (u,v) = iA > (u,v) = (Q\<p{u)<f)i{v)\Q) 

iA—(u,v) = iA\u,v) = (Q\T[<f>{u)<f>\v)]\n) , (2.4) 
and similarly the fermionic Green functions 

is++(u,v) = iS^( Ui v)= (nm a ( u )Mv)]\n) 

iS+p(u,v) = iS<p(u,v) = -(n|$ij(T;)V>a(u)|fi) 

is-p(u,v) = isi (u,v)= (n\T[^ a (u)Mv)}\n) , (2.5) 



where T (T) denotes time ('anti-time') ordering and the additional minus sign in the second fermion 
line is due to the anticommutation property of the fermionic fields. 

These definitions imply immediately the following hermiticity properties for the Wightman functions 

(iA < ' > ( M ,t;)) t = iA^fau) (2.6) 
(i7°S ,< ' > («,v)) t = ^7°5 < - > (w,m). (2.7) 

The four two-point functions (J2.4|) and (J2.5j) are not completely independent. Indeed, with the defini- 
tion of time and anti-time ordering one finds for the chronological (Feynman) and anti-chronological 
Green functions 

G\u,v) = 9(u - v )G > (u,v) + 9(v - u )G < (u,v) 

G\u, v) = 6(u - v )G < (u, v) + 6(v - u )G > (u, v) , (2.8) 

where we have used the generic notation G = {A, S}. From now on we use this notation in relations 
which are identical for bosonic and fermionic Green functions. 

From the definitions of the retarded and advanced propagators 

/^IT /~lt >^Y< >^Y> fit 

Kj = (_t — Kj — \j — Lt 

G a = G* -G> = G< - G l (2.9) 
and the definitions of their hermitean and antihermitean parts 

G h = ^(G r + G a ) 

G a = ^(G a - G r ) = % -{G> -G<) = A, (2.10) 
one obtains the spectral relation 

G h (u, v) = -i sign(u° - v°) A{u, v) , (2.11) 
where sign(w° — v°) = 0(u° — v°) — Q(v° — u°), = (u°, u). A is called the spectral function. 
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2.2. Lagrange density 



We shall now consider the dynamics of fermionic and scalar particles in the presence of a scalar 
condensate which may give space-time dependent masses to both types of particles, and we shall 
assume that scalars and fermions couple via a Yukawa interaction term. The tree-level action is then 

J[0,V] = J d\C (2.12) 

C = # fa - ^ L m^ R - $tfn*i> L + (d^ - f M 2 + £ int , (2.13) 
where the Yukawa interaction Lagrangean reads 

Ant = -^LVHr- 

= -i> (P R ® # + <g> ^ , (2.14) 

and we can rewrite the fermionic mass term as 

-ipLmipR - i'R'm*i'L = — ^(l <S> m ft + ry 5 <g> m a )ij) . (2.15) 

Fermionic fields with definite chirality ipR t i = Pr,l^ are defined with the help of the chirality pro- 
jection operator Pr s l = (1 ±7 5 )/2, and 7 5 = 27°7 1 7 2 7 3 . Our analysis includes the possibility of 
particle flavor mixing, but for notational simplicity we have suppressed scalar and fermionic flavor 
labels (we shall explicitly state so when we consider the special case of the single field dynamics 
without mixing.) Thus, and ip are vectors in the scalar and fermionic flavor space, respectively, 
and y is a vector in the scalar and a matrix in the fermionic flavor space. Then ® denotes a direct 
product of the spinor and fermionic flavor spaces. M and m are complex matrices in the scalar and 
fermionic flavor space, respectively, whose non-diagonal elements couple fields of different flavors to 
each other. Even though the Lagrangean ()2.13)1 cannot fully describe the physics of the Standard 
Model, since it does not contain gauge fields, it contains many essential elements of the physics of 
CP-violation and interactions of importance for baryogenesis. Since the main purpose of this work 
is to establish controlled calculational methods, the Lagrangean (|2.13|) should be a good toy model 
for studying transport, in particular for baryogenesis. 

An important example of the scalar field condensate is the Higgs field at a first order electroweak 
phase transition in variants of the standard model. This phase transition proceeds via nucleation 
and growth of the broken phase bubbles of a nonzero Higgs condensate, which varies at the interface 
of the symmetric and the broken phases, the so called bubble wall. The phase transition bubbles 
quickly grow large in comparison with the thickness of the phase transition front, so for studies of 
local transport phenomena at the scale of the bubble walls one can to a good approximation assume 
planar symmetry. We will explicitly make use of this symmetry in later sections. To allow for 
CP-violation we choose the fermion mass to be a complex matrix 

m(u) = m h (u) +im a (u) = \m(u) |e je(M) , (2.16) 

where rrih and m a denote the hermitean and antihermitean part of m, respectively. The last part of 
this equation is to be understood as an equation for the components. CP-violation can be mediated 
either through m a or complex off-diagonal entries of rrih- The scalar mass matrix is by construction 
hermitean, M 2 ^ = M 2 . Nevertheless, in the case of flavor mixing CP-violation can be mediated 
through the off-diagonal elements of M 2 , provided they are complex. 
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2.3. Two-particle-irreducible (2PI) effective action 



We shall now review how to derive the equation of motion for the Green functions in the CTP- 
formalism. Note first that the tree-level equation can be obtained in a straightforward manner as 
follows. By varying the action (j2.12j) with respect to 0', one obtains the familiar Klein-Gordon 
equation 

(-dl-M 2 (u))<j ) ( U )=0. (2.17) 

After multiplying this from the left by <p'(v) and taking the expectation value one gets the tree- level 
equation of motion for the Wightman function (|2.4j) : 

(-d 2 u - M 2 {u))tA < {u,v) = 0. (2.18) 

By an analogous procedure one finds that the same equation holds for iA > . The equations of motion 
for zA* and zA* are obtained by imposing the appropriate time ordering, such that they acquire on 
the r.h.s. the ^-function sources i5 4 {u — v) and — i5 4 {u — v), respectively. Similarly one finds that at 
tree-level the fermionic Wightman functions (j2.5j) obey the Dirac equation 



(z $ u -m h {u) - rfm a {u)) iS^iu, v) = 0, (2.19) 

and the equations of motion for iS* and iS t again acquire z'5 4 (w — v) and — i5 4 (u — v) on the r.h.s., 
respectively. 

In order to include interactions into the equations of motion, we use the 2PI effective action approach, 
which was originally developed for equilibrium problems in condensed matter physics (where it is 
better known as ^-derivable actions) by Luttinger and Ward (5?} and in relativistic field theory by 
Cornwall, Jackiw and Tomboulis and then adapted to non-equilibrium situations by Chou, Su, 
Hao and Yu |5|, and by Calzetta and Hn 0. In this approach one uses the two-particle irreducible 
(2PI) out-of-equilibrium formulation for the CTP effective action T, which is a functional not only of 
the expectation value of the field <p(u) = (Q\<p(u)\Q) , but also of the two-point functions iA(u, v) and 
iS(u,v) defined in (j2.lH2.2j) . The main advantage of this formalism compared to the more standard 
1PI approach is that, at a given order in the loop expansion, varying the 2PI action probes much more 
accurately the field configuration space. Technically speaking, at a given order in the loop expansion 
the 2PI effective action typically resums a much larger set of diagrams than the corresponding 1PI 
effective action. In absence of sources, the equations of motion are obtained simply by extremising 
the effective action with respect to <f(u), A(u,v) and S(u,v): 



sry,A,s] 

ST[<p,A,S] 

5A 
ST[<p,A,S\ 

5S 






0. (2.20) 



These equations are obtained by the variation of the 2PI effective action, and hence they correspond 
to the Schwinger-Dyson equations, with the self-energy approximated by the 1PI diagrams. Since 
we are here not interested in modeling the dynamics of the scalar field condensate, we shall not 
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study the equation of motion for (p. For our purposes we consider p as given and simply absorb 
the effects of p into space-time dependent mass terms m and M 2 . The quantum dynamics of a 
self-interacting scalar field in the presence of a condensate describing the inflaton decay is considered 
in the framework of the 2PI effective action approach for example in B0ES, while the equivalent 
problem for classical scalar fields is studied in jf^, E3- Various aspects of the out-of-equilibrium 



dynamics and thermalization of quantum fields are considered in [6J, |65|, |6y, |67|. Here we develop 
a gradient approximation for the two-point function, starting with the 2PI effective action, which is 
appropriate for studying the dynamics of fields in the presence of slowly varying backgrounds. 



The 2PI effective action corresponding to the classical action (j2.12j) can be written in the form 

T[A,S] = zTr(A (orl A) - zTr(5 (orl 5) +iTrlnA- 1 - zTr In S' 1 + T 2 [A, S] , (2.21) 

where the minus signs in the fermionic terms are related to Pauli statistics. We use a condensed 
notation, with Tr denoting both integration over space-time and summation over spinor and flavor 
indices. The first two terms in ()2.21j) are the classical (tree-level) actions, the next two terms 
correspond to the one-loop vacuum diagrams shown in figure El while the last term r 2 stands for 
the sum of all two-particle irreducible vacuum graphs. In Paper II we illustrate how to compute 
T 2 in practice. The inverse free propagators A(°) _i and S^' 1 can be read off from the classical 
action (|2.12H2.l'3j) rewritten as 

I[<f> } if>] = f d 4 ud 4 v^(u)A^~\u,v)(f)(v)+ [ d 4 ud 4 vij(u)S {or \u,v)ij(v)+ [ d 4 uC int . (2.22) 



JC JC JC 

They satisfy equations (j2.18M2.lTJ|) and are given by 

A^~\u,v) = (-fig- M 2 (u))5 4 (u-v) 

SW~\u,v) = (z^-m,-*7 5 m a )^( M -^). (2.23) 

We now take the functional derivatives of the effective action ()2.21|) with respect to A and S to 
obtain 

= ^ S <«»- 1 („., t , ) + J S-'( U v,) + §^l = 0. (2.24) 
ob(v,u) oS{v,u) 
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By making use of the definitions for the scalar and fermionic self-energies 

n(«,«) S (2.25) 
oA(v, u) 

. 6T 2 [A,S] 

E(u,v) = -2-^7 r 2.26 

db(v , u) 

and multiplying from the right by A and —S, respectively, we can recast (J2.24)) as the Schwinger- 
Dyson equations 

(-d 2 u - M 2 {u))iA{u,v) = i5 A c {u-v)+ [ d 4 wU(u,w)iA(w,v) (2.27) 



So far we have written everything in the complex contour notation. In the index notation the self- 
energies are 

n°W) = iab 5 ^ S \ (2.29) 

Z«\u,v) = -iab fllf> S] (2.30) 
db Da (v } u) 

and the equations of motion (j2.27H2.2Rjl read 

(-d 2 u - M 2 (u))iA ab (u,v) = aiS ab S 4 (u-v) + J2 c \ d 4 wU ac (u,w)iA cb (w,v) (2.31) 
(ifi u -m h {u) -i^m a (u))iS ab (u,v) = ai5 ab 5 4 (u - v) + / d 4 wY, ac {u,w)iS cb {w,v) .{232) 

c J 

These are the fundamental quantum dynamical equations corresponding to the action (j2.12M2.l3*j) . 
They look deceptively simple, since the full complexity of the problem is hidden in the self-energies , 
which are very complicated functionals of the Green functions and in general not known completely. 

From Eqs. (j2.4M2.5j) and (j2.31M2.32~j) we infer the bosonic equations 

{-d 2 - M 2 )A r ' a - W' a A r ' a = 5 (2.33) 

(-d 2 - m 2 )a<'> - rr © a<'> = n<'> © A a , (2.34) 

and the fermionic equations 

(i0 -m h - rfm a )S r > a - S r ' a S r ' a = 5 (2.35) 
(i$ -m h - z 7 5 m a )5 < ' > - S r 5<'> = £<'> S a , (2.36) 

where stands for integration over the intermediate variable. We use the notation II* = n ++ , etc. 
{cf. Eqs. fl2.4M2.5j) ) and have defined retarded and advanced self-energies in analogy with (j2.9j) . As 
we have already pointed out, not all Green functions are independent: G r and G a can be expressed in 
terms of G < and G > (cf. Eqs. (l2.8H2.9j) ). implying that the equations for the retarded and advanced 
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FIG. 3: Integration contour for the retarded propagators A r and S r in (|2. 38112,3!?)) , 



A a S a 



FIG. 4: Integration contour for the advanced propagators A a and S a in (|2.38H2.39|) . 

propagators are redundant. It is not difficult to show that Eqs. (j2.33H2.3B|) are consistent provided 
the relations 

lf(u,v) = 5 4 (u-v)Tl ss (u) + 6(u -v )Tl > (u,v) + 6(vo-uo)Tl < (u,v) 

U\u,v) = 5\u - v)Tl ss (u) + 6(u - vo^fav) + 6(v - u^lPfav) , (2.37) 

as well as corresponding relations for the fermionic self-energies, are satisfied. This should be the 
case for any reasonable approximation of the self-energies. The additional singular terms I1 ss (m) and 
£ ss (m) we have allowed for appear naturally in some theories and can be absorbed by the mass terms 
(mass renormalization) . An example of a singular self-energy contribution is the one-loop tadpole of 
a scalar theory with quartic self-interaction. The one-loop expressions for the self-energies that we 
will derive in Paper II indeed satisfy ()2.37|) and the corresponding fermionic equations. 

To clarify the physical meaning of Eqs. (j2.33M2.3^j) . we observe that the equations for the retarded 
and advanced propagators ()2.33j) and ()2.35j) describe mostly the spectral properties of the system. 
Indeed, this can be seen for example from the tree-level solutions, which in the limit of constant mass 
terms and no flavor mixing read 

A °>^ = J ^ HU - \^M^L M e (2 ' 38) 
S r > a (uv)- I d ' k r-^-(u-v) ^ + m h -z^m a 

b ° {%V) ~ J (2vr) 4e k 2 -ml-ml±tega(ko)e' 1 J 

where we use the infinitesimal pole shifts =Fze (e — > 0+) to indicate the standard integration prescrip- 
tion for the retarded and advanced propagators, respectively, shown in figures H3 and 0] The spectral 
functions (j2.1U|) are then 

A* = j 0^e- ik < u -^sign(k o )6(k 2 - M 2 ) (2.40) 

-V = j -0^e- iHu - v \ft + m h - i^m a )^ign{h)5{k 2 -m\- m 2 a ) . (2.41) 

The (^-functions indicate that the frequencies of the plasma excitations are constrained to lie on the 
mass shell, which for scalars is given by 



±w , u^=\lk 2 + M 2 (2.42) 
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and for fermions by 

k = ±u , lu = Jp + m 2 h + ml . (2.43) 



On the other hand, the equations of motion for the Wightman functions ()2.34|) and (j2.36|) describe 
mostly statistical (kinetic) properties of the system. This can be seen for example from considering 
the thermal equilibrium solutions to these equations, which are given in section l2~5l below. The kinetic 
equations ()2.34j) and ()2.3()j) can be rewritten in the form of the Kadanoff-Baym (KB) equations 

{-d 2 - m 2 )a<'> - n h © a<'> - n<'> © A ft = c^^(n > 0A < -n < A>) (2.44) 

(^ -m h - i 7 5 m a )S ,< ' > - Eft - £<'> S h = = ^(S> S< - S< S>) , (2.45) 

where and denote the scalar and fermionic collision terms, respectively, and the hermitean 
parts of the self-energies are defined analogously to those for the Wightman functions ()2.10|) . With 
reasonable approximations for the self-energies, the Kadanoff-Baym equations are suited to study ki- 
netics of out-of-equilibrium quantum systems. The terms IIftA < and Sft5' < on the l.h.s. of f)2.44H2.4H|) 
represent the self-energy contributions, which are sometimes considered as a (nonlocal) contribution 
to the mass terms. We postpone the discussion of how inclusion of the self-energy may change the 
plasma dynamics and in particular the quasiparticle picture to a future work. The collision terms 
and C,p, as we will see in Paper II, contain the standard gain and loss terms responsible for equilibra- 
tion, but they may also contain additional CP-violating sources. Finally, the terms II < Aft and S < 5'ft 
induce broadening of the on-shell dispersion relations and may cause breakdown of the quasiparticle 
picture. We shall consider the role of these terms in some detail in sectional There we will need the 
width 

T^ = -(Tl a -W)= l -(Tl>-n<) (2.46) 
and the hermitean part of the scalar self-energy 

lift = ^(ir + it) = iT - l(n> + n<) , (2.47) 

which satisfy the spectral relation (cf. Eq. (j2.11D ) 

U h (u, v) = ^sign(u° - v°) [IP (it, v) - n<(u, v)} = -isign(w° - v°) T^u, v) . (2.48) 
Analogous relations hold for the fermionic self-energies. 

2.4. Wigner representation and gradient expansion 

In equilibrium the Green functions A(u,v) and S(u,v) depend only on the relative (microscopic) 
coordinate r = u — v. This dependence corresponds to the internal fluctuations that typically take 
place on microscopic scales. In non-equilibrium situations, however, A and S depend also on the 
average (macroscopic) coordinate x = (u + v)/2. This dependence describes the system's behavior 
on large, macroscopic scales. For example, if the system couples to an external field which varies in 
space and time, the Green functions will show a corresponding dependence on the average coordinate. 
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Thick bubble walls at a first order electroweak phase transition represent precisely an example of 
such an external field. 

In order to separate the microscopic scale fluctuations from the behavior on macroscopical scales 
one typically performs a Wigner transform, which is a Fourier transform with respect to the relative 
coordinate r. The Green function in the Wigner representation, called Wigner function, is 

G(k, x) = J d 4 r e ihr G(x + r/2, x - r/2) . (2.49) 

In the Wigner representation the hermiticity properties ()2.6|) and (J2.7)) become simply 

(tA < ' > {k,x)) j = iA^fcx) (2.50) 
(rfS^fcx))* = tfS^fcx). (2.51) 

In order to transform the equations of motion, we make use of the general relation 

d\u - v) e ik < u ~ v ) J d 4 wA(u, w)B(w, v) = e~ i0 {A(k, x)}{B(k, x)} , (2.52) 

where x is the average coordinate and the diamond operator is defined by 

o {■}{■} = \{d {1) ■ d ( * } - ■ d^) {•}{•}. (2.53) 

The superscripts (1) and (2) refer to the first and second argument, respectively, and d = d x . Since 
the nonlocal terms have the form of 1)2.52)1 . it is a simple exercise to show that the Kadanoff-Baym 
equations (j2.44H2.45*)) . when written in the Wigner representation, become 

(k 2 -^d 2 + tk-d- M 2 e"^) A<'> - e- io {n4{A<' > } - e-^II^HAfc} = (2.54) 

-mne-fr 8 * - i^m^^S^ - e^{E ?l }{^ < ' > } - e^{^>}{S h } = , (2.55) 
where the collision terms are 

C = ie-*({n>}{A<} - {n<}{A>}) (2.56) 

Cip = ie"*({E>}{5<} - {S<}{5>}) . (2.57) 

We shall analyze these equations in gradient expansion: we assume that the variation of the back- 
ground field, and therefore also the variation of the Green functions and self-energies, with respect 
to the macroscopic coordinate x is small when compared with the momenta of plasma excitations. 
Then we can perform an expansion in derivatives with respect to x. Formally, this criterion reads 

d < k. (2.58) 

That is, the background field is assumed to have a characteristic length scale that is large in com- 
parison to the de Broglie wavelength of the particles in the plasma. 



This assumption of a slowly varying background field seems to be justified in our case [681] : in the 
MSSM, for instance, the width of the bubble wall L w is roughly 10/T, where T is the temperature 



14 



of the plasma. The typical momentum of a particle in the plasma is of the order T, so that the 
de Broglie wavelength I^b ~ 1/T is indeed small when compared to L w . Since the expansion in 
powers of derivatives can be viewed as an expansion in powers of the Planck constant h, the gradient 
expansion is equivalent to a semiclassical expansion. We expect that the leading order terms in the 
gradient expansion correspond to classical kinetic equations, while higher order derivatives represent 
quantum corrections. 



2.5. Sum rules and Equilibrium Green functions 



It is easily shown that, as a consequence of the equal-time commutation relation for scalars, the 
Wigner transform of their spectral function (j2.10|) obeys the spectral sum rule 

/oo JI, 
^k Q A 4> (k,x) = l, (2.59) 
■oo K 



while the anticommutation rule for fermions leads to 

-'fA i> (k,x) = l (2.60) 



T- 

for the fermionic spectral function. These sum rules are important, since they have to be imposed 
as consistency conditions on the solutions for the Wigner functions G K and G > , or for G a and G r , 
as indicated in (j2.10j) . Furthermore, making use of these spectral sum rules and the Kubo-Martin- 
Schwinger (KMS) relations 

A>(&) = e^A<(A;) 

5>(*) = -eF*S<(k), (2.61) 

where (3 = 1/T denotes the inverse temperature, one can unambiguously obtain the thermal equi- 
librium propagators 0|. They are translationally invariant in space and time and hence, when 
written in the Fourier space, defined by the Fourier transform with respect to the relative coordinate 
r = u — v, they display a dependence on the momentum only: 

G eq (k) = J d*re tk r G cq (r). (2.62) 

The bosonic and fermionic equilibrium propagators are the solutions of the Klein-Gordon 1)2.18)1 and 
the Dirac equation (J2.19)) . respectively, which in the Fourier space read 

(k 2 - M 2 )zA< q >(A;) = , (k 2 - M 2 )iA r £{k) = i 

(f-m h - % ^m a ) tS< q >(k) = , {f-m h - i^m,) tS r eq a (k) = i , (2.63) 

with the normalization (j2.59H2.6U|) and the boundary conditions (|2.61|) . The scalar equilibrium Green 
functions are 

iAl q (k) = fc2 _ M2+ \ sign(fco)6 + 2n5(k 2 - M 2 )sign(k )nt q (k ) (2.64) 

^Al q (k) = - fc2 _ M2+ ' jsign(fcb)e + 2 ^( fc2 " M 2 )sign(fc )(l + <(^o)) (2.65) 

iA<(k) = 2n5(k 2 - M 2 )sign(k )ni q (k Q ) (2.66) 
iA> (k) = 27c5(k 2 -M 2 ) S ign(ko)(l+ni q (k )), (2.67) 
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where the Bose-Einstein distribution appears: 



<(**) = ^731 • ( 2 - 68 ) 



Similarly, for the fermions we have 



iS^k) = ^7T1 - MP + m h~ rfm a )5(k 2 - \m\ 2 )sign(k )n cq (k ) (2.69) 

/C \7Tl\ ~r~ 2S1^I1( /Cq )€ 

^<(A;) = -2tt(^ + m h - i 1 5 m a )5(k 2 - |m| 2 )sign(A;o)n cq (A; ) (2.71) 

iS>(k) = 2ir0 + m h -i 1 5 m a )5(k 2 -\m\ 2 )sign(k o )(l-n cq (k o )), (2.72) 

with the Fermi-Dirac distribution function 

and \m\ 2 = m 2 + For the case of several mixing particle species, the equilibrium Green functions 
are diagonal matrices in flavor space, and the given relations hold for the diagonal elements. 
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3. REDUCTION TO THE ON-SHELL LIMIT 



In this section [79] we pay particular attention to the self-consistency of the on-shell limit as a 
controlled expansion in the coupling constants of the theory. To this aim we consider the role 
of the dissipative terms on the l.h.s. of the propagator equations ()2.33|) and the Kadanoff-Baym 
equation (|2.54p . We shall see that including these terms into the retarded and advanced propagators 
results in a Breit-Wigner form for the spectral function. By making use of the spectral decomposition 
ansatz for the Wigner function with the Breit-Wigner form for the spectral function, we shall study 
under which conditions the kinetic equation reduces to the on-shell form. 

For simplicity, we consider here the scalar case. It is quite plausible that similar techniques can be 
used to extend the analysis presented in this section to the case of fermions. Let us begin by rewriting 
equations ()2.33|) and (|2.54j) more compactly, 

e-*{n}±iIV}{A*"} = 1 (3.1) 

e-*{n}}{A<'>} - e-^n^HA*} = , (3.2) 
where C$ is defined in Eq. (I2.56j) . and we defined 

^1 = k 2 -M 2 -U h . (3.3) 

Upon adding and subtracting the complex conjugate of (|3.1|) . we arrive at the following propagator 
equations for the case of a single scalar field up to second order in gradients 

coso{fi}±ir4{A r ' a } = 1 (3.4) 
sino{fi}±ir4{A r * a } = 0, (3.5) 

where we used the fact that A' r,a commute with and T^. This commutation property is special to 
the single field case however. Indeed, in the case of mixing scalars, A r,a are matrices which in general 
do not commute with fi? and T^. 

Now observe that Eq. ()3.5|) can be obtained from Eq. ()3.4|) by an application of the differential 
operator tan o, and hence it gives no new information. Because cos o is an even function, we then see 
that the corrections to the propagator equation appear first time at the second order in gradients |7o| . 
Therefore, to the first order in gradients, we have simply 

A^ = (^±irV)-\ (3.6) 

or equivalently (c/. (|2.10|) ) 

^ = 7^72 (3-7) 

iL <t> i" L </> 

The kinetic and constraint equations can be obtained from the Kadanoff-Baym equations ()3.2|) as 
the hermitean and antihermitean parts, 

-o{nJ}{iA<'>}+o{ffl<->}{A fc } = ^(n>A<-n<A>). (3.9) 
^A<'> + %({n>}{A<}-{n<}{A>}) = n<A„ (3.10) 
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where we truncated to second and first order in gradients, respectively. Because of the hermiticity 
property (|2.5()jl of the scalar Wigner function, both kinetic and constraint equations must be satisfied 
simultaneously The constraint equation selects the physical set of solutions among all solutions of 
the kinetic equation. 

Note that the equations (j3.7H3.8|) are only formal solutions in the sense that and are func- 
tional of iA <,> and iS <y> . Consequently, equations (J3.9H3.10|) are complicated integro-differential 
equations. However, in many cases deviations from equilibrium are small, so that one can linearize 
them in deviation from equilibrium, such that Ut and can be computed by using the equilibrium 
distribution functions (|2 . 64112. 73|) . 

Inspired by the spectral form of the equilibrium solutions, we shall assume that close to equilibrium 
the following spectral decomposition holds 

iA < = 2A 4 ,n^ tA> = 2^(1 + (3.11) 

where = n^{k : x) is some unknown function representing a generalized distribution function on 
phase space. Similarly, for the self-energies, we can formally write 

zn< = 2r>n, ffl> = 21^(1 + n n ). (3.12) 

We have adopted a somewhat nonstandard definition for Tf, the more standard definition is obtained 
by the simple replacement — > 2k T l f ) . At the moment % is just an arbitrary variable replacing 
n <,> , but it will become closely related to the particle distribution function in the end. It is important 
to note that unlike n^, it is not to be considered a free variable in the equations. 
Then using the definitions (j3.11H3.r2*j) and the identity 

o{f}{gh} = go {f}{h} + ho {f}{g}, (3.13) 

we rewrite Eqs. f|3.9H3.10|) as 

A^o{^}{n4 + T (/> o{A h }{n^} = T^A^ - n u ) (3.14) 
ttlA^nj, + o{T, i) A ( j > }{n ( j ) } = T^riuAh. (3.15) 

We can now see that taking nn — > solves the kinetic equation 1)3.14)1 to the zeroth order in 
gradients, establishing the promised connection between nn and n^. This also shows that, to the 
first order in gradients, it is consistent to replace nn by the zeroth order solution n^, whenever the o 
operator acts on nn- In fact, we already took this into account by replacing nn by in the T^-terms 
on the l.h.s. of Eqs. (J3.14M3.l5|) . 

We shall now see how to reduce these equations to an on-shell approximation, and in particular ask to 
what extent the on-shell limit is a good approximation to the plasma dynamics close to equilibrium. 
The equations can be simplified further when the width is small. While this is in general true in 
the weak coupling limit, due care must be exercised in making this approximation, because both 11^ 
and are often controlled by the same coupling constants. 
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and to first order in gradients, the spectral function (J3.8)) reduces 



3.1. Propagator equation 

In the weak coupling limit 
to the singular spectral form, 



A4, >• A s = vr sign(/c ) 6(pfy, 



(3.16) 



where the correction is of order T^. This singular hypersurface defines the usual quasiparticle dis- 
persion relation 

Ql = k\ -k 2 -M 2 - n^fco, k, x) = 0, (3.17) 

which defines the spectrum of the quasiparticle excitations of the system. This equation has in 
general two distinct solutions ko = u>k and ko = —Uk, corresponding to particles and antiparticles, 
respectively. Thus the spectral function (jH.lfij) breaks up into two clearly separate contributions: 

7T 7T - 

A s = - — Z k 5(k - LO k ) - ^—Zk 5(k + u k ), (3.18) 

ZUJk liO k 

where the wave function renormalization factors Z k and Z k are 

2u k Z^ = \dnl/dk Q \ h)=uk , 2u k Z^ = \dnl/dk \ ko= _ Qk . (3.19) 

The definition of in (|3.7|) becomes more problematic in the limit 1^ — > because of the appearance 
of the so called Landau ghost poles, when perturbative expressions are used for the self-energy LT^ [70|. 
Therefore, when needed, we shall relate to A<j, via equations (j3.4H3.5|) before taking the zero width 
limit. Let us note however that the spectral representation 

A h (k,x) Re f dk' ^lM. (3.20) 
J k -k' + te 

remains valid for A^ even in the on-shell limit [to| . 

The full spectral function A4 must obey the sum- rule (|2.59|) . The fact that the singular function 
A s ()3.18|) fails this rule when interactions are included, tells us something important about the 
physical nature of the quasiparticle approximation. Indeed, one is here describing the spectrum of 
inherently collective plasma excitations by a simple single particle picture. The amount by which 
A s fails the sum-rule can then be attributed to the effect of collective plasma excitations. This 
gives a quantitative estimate for the goodness, but also indicates the limits of applicability of the 
quasiparticle approximation. 

Since both 11^ and are controlled by the same coupling constants, they typically acquire contri- 
butions at the same order in the coupling constant, so that one cannot take the limit — > as an 
independent approximation. The theories containing tadpoles, such as the A0 4 -theory, seem at a first 
sight to be an exception. However, the tadpoles give rise to a singular non-dissipative contribution to 
the self-energy, which can be absorbed by the mass term, implying that they cannot be responsible 
for thermalization. This means that in order to study thermalization properly, one has to include 
the higher order nonlocal dissipative contributions, which typically contribute to both the self-energy 
and width at the same order in the coupling constant. Consequently, there is a danger that one 
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might lose the consistency of the single particle picture. Nevertheless, as we shall now see, this is not 
the case. Namely, the two quantities appear in the spectral function in quite distinct ways such that 
the kinetic equation allows a consistent on-shell limit assuming equal accuracy in the computation 
of Il/j and r^,. 



3.2. Kinetic equation in the on-shell limit 



Our aim is now to study the conditions under which the kinetic and constraint equations (J3.14I - 
13.15)1 reduce to the on-shell approximation for the effective plasma degrees of freedom found in the 
previous section. Bearing in mind the discussion concerning consistency of the quasiparticle limit, 
we start by using the expressions (|3.7H3.8|) with a finite width T^. The term linear in on the l.h.s. 
of (j3.14|) is often simply dropped without a careful consideration £lj|. However, making use of the 
identity T^A^ = f&A^,, which is accurate to the leading order in gradients (cf. Eq. ()3.15|) ). and 
expressions (J3.7H3.8j) . we can in fact combine it with the first term in the kinetic equation to obtain 

— 2r 3 2T 2 Q 2 
( fl 4 + ]4)2 ° WIN) - (Q4 I r t) 2 ° ( r 4 W) = rV^K, - n n ) . (3.21) 

Note that in the weak coupling limit and to leading order in gradients 



-2Tl 



° A s , (3.22) 



such that Eq. (|3.21|) can be rewritten as 

A s o {«}}{n*} - 2T (t> A s A h o {T^}{n^} = T^A^ - n n ) . (3.23) 

Under the reasonable assumption that is a smooth function of the energy, we get from ()3.23|) 
a consistent expansion of the final result in powers of evaluated at the pole. In this expansion 
the contributions from the off-shellness cancel so that up to order 0(T^) the integrated equation 
involves only the leading on-shell excitations, defined by the singular spectral function ()3.18j) . In 
other words, the integrated kinetic equation exhibits the important property of closure onto the 
on-shell excitations in the weak coupling limit. We can then replace ()3.23|) by the singular on-shell 
kinetic equation 

A«{0;}{n4 = r^(n r n n ). (3.24) 

We again emphasize that it is actually consistent to include non-local loop contributions to 11^ on 
the l.h.s. despite the fact that the singular single particle limit does not, strictly speaking, exist for 
the excitations. This is simply because the corrections due to off-shellness cancel up to 0(T$). This 
is fortunate, because applying the weak coupling limit in the strict sense to (I3.24|) indeed requires 
either computing TE^ and to the same order, or simply neglecting T^, leading to a collisionless 
Boltzmann equation. A somewhat orthogonal approach has been taken by Leupold [7^ in which, 
based on particle number conservation, the author advocates a modification in the relation between 
the Wigner function and the on-shell distribution function. 

The final step is to show that the constraint equation (j3.15|) can be written as 

njAK - n n) = - o {IV4 S }{™4- ( 3 - 25 ) 
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In the limit 1^ — > the constraint equation becomes identically solved when the quasiparticle on- 
shell condition (j3.17j) is satisfied. In other words, the on-shell condition and constraint equations 
become degenerate. Note that we can here refer to the on-shell condition, not necessarily restricted 
to the singular approximation, but in the broader sense we derived the on-shell kinetic equation ([3.24)1 . 
Indeed, the constraint equation can always be satisfied by the on-shell solutions as given by Eq. (|3.24|) . 
whereby it only gives rise to an integral constraint for the off-shell excitations. These however, do not 
belong to the set of dynamical degrees of freedom for the kinetic equation (|3.24|) . which only contains 
the on-shell excitations as defined by the singular overall projection operator A s . In this sense the 
kinetic equation (|3.24J) and the constraint equation decouple. Of course, the off-shell excitations in 
the constraint equation ()3.25|) become relevant when the off-shell effects are included into the kinetic 
equation. 

We emphasise that the derivation of the on shell kinetic equation (|3.24|) applies to the special case 
of one scalar field in weak coupling limit and close to equilibrium. It would be of interest to extend 
the analysis to the case of mixing scalar and fermionic fields. 
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4. KINETICS OF SCALARS: TREE-LEVEL ANALYSIS 



In this section we analyze the tree-level dynamics of the scalar sector of our theory with mixing and 
emphasise its ramifications. The hermiticity property ()2.50|) of the scalar Wigner function implies 
that the hermitean and the antihermitean part of the equation of motion 



k 2 --d 2 + ik-d 
4 



M 2 (x)e-^ d ' dk 







(4.1) 



ought to be satisfied simultaneously. A simple analysis reveals that the hermitean part corresponds to 
the kinetic equation, while the antihermitean part corresponds to the constraint equation. Roughly 
speaking, the kinetic equation describes the dynamics of quantum fields, while the constraint equation 



constrains the space of solutions of the kinetic equation 52]. In the simple case when N — 1 it is 
immediately clear that the first quantum correction to the constraint equation, the real part of ()4.1|) . 
is of second order and to the kinetic equation (imaginary part) of third order in gradients or second 
order in h. To extract the spectral information to second order in h is quite delicate since the 
constraint equation in ()4.1|) contains derivatives. The situation is more involved in the case of more 
than one mixing field. In this case it is convenient to rotate into the mass eigenbasis 



M 2 d = UM 2 U\ 



(4.2) 



where U is the unitary matrix that diagonalizes M 2 . In this propagating basis equation (|4.1|) becomes 



(k 2 -h) 2 + ik-V- M 2 e-^-^j A< = 0, 
where A^ = UA < U^ and the 'covariant' derivative is defined as 



V ll = d^-i [E^ 



E u = iUd„U ] . 



(4.3) 



(4.4) 



i^V, = 2VA<, and M 2 cJ A< = 



Making use of (iA|) f = iA^, = which implies (V^iA^Y 
^{M 2 S , iAj} + ^[M 2 s , iA^], we can identify the antihermitean part of (|4.3|) as the constraint equation 



k 2 -\v 2 yA< -\{MlrA< 



M' zA< 



0. 



and the hermitean part is the kinetic equation 



k.V l A<+ 1 -{M 2 , i A<}+'- 



M?,iA' 



0. 



(4.5) 



(4.6) 



We defined 



M 2 C = M 2 cos - V ■ d k 
M 2 = M 2 d sin - V ■ d k . 



(4.7) 



and [-,-], { ■ , ■} denote a commutator and an anticommutator, respectively. The constraint (|4.5j) 
and kinetic equation ()4.6j) are formally still exact. Since we are interested in order h correction to 
the classical approximation, it is good to keep in mind that one can always restore h dependencies 
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by the simple replacements d — > hd and iA^ — > H~ l iA^. We work to order h accuracy, so it suffices 
to truncate the constraint equation (|4.5|) to first order in gradients, and the kinetic equation (|4.fi|) to 
second order in gradients: 

kHA< - \ {M 2 d ,tA<} + I [DM d 2 ,4«A<] = 0, (4.8) 



fc-P^A< + i{PM d 2 ,^A<} + ^ 



M](l-i(B-a fe ) 2 ),iA< 



0. (4.9) 



We shall now consider these equations in more detail. The off-diagonal elements of iA^ in both 
the constraint ()4.8j) and the kinetic equation ()4.9j) are sourced by the diagonal elements through 
the terms involving commutators, which are suppressed by at least h with respect to the diagonal 
elements 0. On the other hand/the off-diagonals source the diagonal equations through terms that 
are of the same order as the diagonals. For the CP-violating diagonal densities that are suppressed 
at least by one power of H, this immediately implies that the CP-violating off-diagonals are at least 
of the order ti 2 , and thus cannot source the diagonal densities at order h. By a similar argument, 
the second order term in the commutator in ()4.9|) can be dropped, since it can induce effects only at 
order h 2 . 

This analysis is however incomplete for the following reason. When the rotation matrices = 
iUdfjJP in (|4.4j) contain CP-violation, the off-diagonal elements can in principle be CP-violating 
already at order h, and hence source the diagonals at the same order. In order to get more insight 
into the role of the off-diagonals, we now analyze the case of two mixing scalars. For notational 
simplicity we omit the index d for diagonal in the following. All quantities have to be taken in 
the rotated basis, where the mass is diagonal. The constraint equations ()4.8|) . when written in 
components, are 

(V - M^iA< =f \5(M 2 ) (~ 12 • d k iA< + S 21 • d k iA< ) = 0, (4.10) 

k 2 -M 2 + l -(d5(M 2 ))-d k yAf 2 + h(M 2 )Z 12 -d k 5(iA<) = 0, (4.11) 

where M 2 = Tr(M 2 )/2 = (Mf 1 + M| 2 )/2, 5(M 2 ) = M\ x - M 2 2 , 5(iA<) = zA< - iA< 2 , and the 
equation for iA 21 = (iAf 2 )* is obtained by taking the complex conjugate of (|4.11|) . To order h the 
diagonal equations ()4.1U|) are solved by the spectral on-shell solution 



iA£(Jfe,z) = 27r(5^ 2 -M 2 1 )sign(r€o)^f(A;,x) 

= —[5(k -uj l f ) i)-8(ko + u < j ) i)]nt(k,x), (4-12) 

where = [k 2 + M^] 1 / 2 , and nf(k,x) represents the occupation number density on phase space, 
which in thermal equilibrium reduces to nf(k,x) — > r^ q = l/(e^ fc ° — 1). By making use of the sum 
rule (|2.59|) for the spectral function A$ = (i/2)(A > — A < ), one can show that the other Wigner 
function has to be of the form 

iA>(fc,x) = 2n5(k 2 - M 2 1 )sign(A;o)(l + nf (k, x)) (4.13) 



with the same density nf as in ([4. 12)1 . These solutions are consistent provided the off-diagonals iAf 2 
and iA 21 are of order h, which, as we shall argue, is a self-consistent assumption. The off-diagonal 
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constraint equation (j4.11j) is solved by 



zA< = n* 5(A; 2 - M 2 ) - ~H 12 • «9 fc Tr(*A<) - fc • S ia 5(zA<) , (4.14) 

where n*f 2 = nf 2 (k,x) is a function that can be determined from the boundary conditions. Since we 
are interested in situations close to equilibrium in which the off-diagonals are driven away from zero 
primarily by the diagonals, we can set nf 2 to zero for our purposes. Note that the solution (|4.14jl 
contains derivatives of the delta function, and hence it does not strictly speaking represent an on-shell 
form. 

Consider now the kinetic equations (|4.9jl . which in components read 

±i k ■ fs 12 iA< - E 21 iAf 2 J 

l -5(M 2 ) (e 12 • d k iA< - E 21 • d k iA< ) (4.15) 

-ifc ■ E 12 5(iA<)-U{M 2 )E 12 ■ d k Ti{iA<) , 

(4.16) 

and the equation for iA 21 is again obtained by taking the complex conjugate of ([4.16)1 . To leading 
order in gradients the off-diagonal equation (|4.16J) is solved by 

<A< = -\ S 12 • «9 fe Tr(*A<) - k ■ H 12 5(zA<) , (4.17) 

and similarly for iA 21 . Remarkably, this corresponds precisely to the constraint equation solu- 
tion (|4.14jl . provided one sets nf 2 = 0, representing a very nontrivial consistency check of the kinetic 
theory for mixing particles. By making use of the technique of Green functions, in Appendix [X] we 
show that the leading order result ()4.17j) represents a valid approximation, provided the condition 
for the gradient approximation k ■ d <C 5(M 2 ) is satisfied, that is one is not near the degenerate mass 
limit, 5(M 2 ) = 0. 

The example of mixing scalars represents indeed a nice illustration of the workings of the constraint 
and kinetic equations. Unlike in the one field case, in which the kinetic flow term can be obtained 
by acting with the bilinear o operator on the constraint equation, in the mixing case the situation is 
more complex. All of the solutions of the kinetic equation 1)4.16)1 are simultaneously solutions of the 
constraint equation (|4.11|) . The converse is however not true: the constraint equation contains a larger 
set of solutions, which include the homogeneous solutions (|4.14|) that lie on a different energy shell, 
which are in fact excluded by imposing the kinetic equation ()4.16j) . This can be shown as follows: 
the constraint equation ()4.11)1 can be obtained from the kinetic equation 1)4.16)1 by multiplying it by 
{k 2 — Tr(M 2 )/2)/5(M 2 ), followed by a partial integration. The constraint equation obtained this 
way has of course a larger set of solutions. The additional solutions are of the type: a function 
multiplying 5(k 2 - Tr(M 2 )/2). 

Finally, upon inserting ()4.17jl into the diagonal equation ()4.15j) . we get 

(k-d+ \{dM 2 ,) -d k ±tk- [S 12 , 3 21 ] • d k )zA< = . (4.18) 



(* • d + -{dMD ■ d k )iA< = 
(k. d+ l -{dM 2 ) ■ d k + U{M 2 )-ik ■ <?(S))iA< = 
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It is remarkable that not only the off-diagonals have disappeared from (|4.18|) . but also the diagonals 



decouple. The commutator k 
and in the wall frame, 



— 12? — 21 



• dk does not vanish in general. However, for planar walls 



k ■ 



— 12? —21 



■d. 



k?. I — Z 12 —z21 — — *21 — 



z21 — zl2 



d k =0 



vanishes identically, so that Eq. (j4.18|) reduces to 



fc.a-i(s,jwS L )^,)/A il = o. 



(4.19) 



(4.20) 



This proves that, although a naive analysis of Eqs. (j4.8H4.9|) indicated that there might be a CP- 
violating source in the kinetic equation at order h, a more detailed look at the structure of these 
equations shows that for planar walls considered in the wall frame no such source appears. Equa- 
tion ()4.20|) thus represents a self-consistent kinetic description of scalar fields accurate to order h, 
such that the only source in the kinetic equation is the classical force, Fi = —Vuj^ = — VM^/2a;^. 
In addition, we have shown that to order h the on-shell approximation for the diagonal occupation 
numbers (14.12)) with the classical dispersion relation fc — i^^j still holds. Even though the on-shell 
approximation fails for the off-diagonals, this is only relevant for the dynamics at higher orders in 
gradient expansion. Further, inclusion of collisions in the off-diagonal equations cannot change our 
conclusions concerning the source cancellation expressed in Eq. ()4.18|) . even though it may introduce 
new collisional contributions. (For a detailed study of collisional sources we refer to Paper II.) Need- 
less to say, this analysis easily generalizes to the case of N mixing scalars. This completes the proof 
that was originally presented in Ref. [^, which states that in the flow term of scalars there is no 
source for baryogenesis at order h in gradient expansion. 



4.1. Boltzmann transport equation for CP- violating scalar densities 

For completeness, and to make a connection between the tree-level analysis presented here and the 
analysis of the collision term in Paper II, we now show how - starting with the kinetic equation 
for the Wightman function ()4.2())1 - one obtains a Boltzmann equation for the CP-violating scalar 
particle densities. All flavor matrices are to be taken in the basis where the mass is diagonal. 
In the beginning we review some of the basics of the C (charge) and CP (charge and parity) trans- 
formations of quantum scalar fields. Under C and P, a scalar field <fi transforms as 

<j?(u) = V <j)(u) V" 1 = rj^u) , V <p\u) V~ l = rj^iu) , (4.21) 

where |^| = 1 and \r]^\ = 1 are global phases, and = (u , —u) denotes the inversion of the spatial 
part of u. Note that our definition of charge conjugation includes an additional transposition with 
respect to the usual definition |73], which is required in the case of mixing scalar fields, when is a 
vector in flavor space. From ()4.21|) it follows that the scalar Wightman functions transform as 

iA < (u,v) -fU iA > (v,u) T (4.22) 
iA<(u,v) ^ tA > (v,u) T . (4.23) 
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When written in the Wigner representation, the equivalent transformations are 

iA < (k,x) iA > (-k,x) T = iA c< {k,x) (4.24) 
tA < {k,x) m> iA > {-k,x) T = iA cp< {k,x) . (4.25) 

The CP-conjugate of the Wigner function iA < (k, x) is related to physical objects at position — x with 
momentum —k. Therefore we introduced additional inversions of the spatial parts of the position and 
momentum arguments in the definition of iA cp< (k,x), so that this object indeed describes particles 
at position x and with momentum k. Note that this definition differs from what is usually found 
in textbooks. In order to study the effects of the C and CP transformations on the distribution 
functions, we proceed analogously to ()4.12|) and define the spectral solutions 

i At p< (k, x) = — [8(k - u+i) - 8(k + uj^)] nf c/cp (k, x) (4.26) 

iAf p> (k, x) = — [5(k - un) - 5(k + u^)] (l + nf c/cp (k, x)) . (4.27) 

Now upon inserting this and (j4.12H4.r3j) into (j4.24j) and (j4.25j) . we find 

n\ c {tjOfr, k,x) = nf cp (lu^, k,x) = - [l + nf ( - u^, -k, x)] . (4.28) 

Because of the additional inversions in the definition of iA cp< , the two densities are identical. 

We shall now consider the kinetic equation for scalars. First we include the collision term into 

Eq. (jonf . 

(k-d- 1 - {d z Ml)d kz ) iA< (*, x) = \ [C#i (k,x)+ Cj ±i (k, x)] , (4.29) 

where Cpa denotes the diagonal elements of the collision term defined in Eq. (j2.56j) after rotation 
into the basis where the mass is diagonal. Next, we insert the spectral on-shell solution (j4.12jl and 
integrate over positive frequencies, to obtain 

d t + — - {dz ^ iiz)) d kz ) ff + (k,x) = ±- [°°dk [C (f>ii (k,x)+Cl ii (k,x)] , (4.30) 
io^i Zuj^i J 2n J 

which is the kinetic equation for the scalar distribution function 

ff + (k,x) =nf(uj^(k } x) } k,x) . (4.31) 

Since iA > (k, x) and iA < (k,x) satisfy identical equations, we can immediately write 

(-&-a+i(3 z ;i4(z))^ (4.32) 

where we took account of (j4.25j) . We can omit the transposition in the collision term, because we are 
only interested in the diagonal elements. Upon integrating over positive frequencies this becomes 

% + — - i9z ^ Mz)) d kz ) ft(k,x) = 1- [°°dk [ C(Pii (-k,x)+Cl ii (-k,x)] . (4.33) 
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This transport equation is the CP-conjugate of (|4.3U|) for antiparticles, where we defined the antipar- 
ticle distribution function as 

ff_(k, x) = nf cp {u^ k,x) = - [1 + nf ( - (J*, -k, x)] . (4.34) 

This then implies the following definition for the CP-violating particle density: 

5f? = ft - ft ■ (4.35) 
The relevant kinetic equation for 5ff is obtained simply by summing (J4.3Uj) and (|4.33|) . 

(% + — - {dz ^ Uz)) d k \ 5ff(k,x) = J- [ O °dk [C (f>ii (k,x)+Cl ii (k 1 x) 

+e rt>ii (-k,x)+el ii (-k,x)]. (4.36) 

It should be stressed that the absence of a CP-violating force in the flow term at linear order in 
h is made explicit in this equation. Further, this equation makes it transparent how to extract 
CP-violating contributions from the collision term. 



4.2. Applications to the stop sector of the MSSM 



This analysis is relevant for exam ple for a calculation of the CP- violating force in the stop sector 



q= (t Ll t R ) T of the MSSM 
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in which the mass matrix reads 



m Q 

y(AH 2 + f M*H 1 ) 



y(A*H 2 + f xH 1 ) 



my 



(4.37) 



where tuq and mfj denote the sum of the soft SUSY-breaking masses, including D-terms and 



y 2 H 2 . Our analysis immediately implies that, for squarks in the quasiparticle picture, there is no 
CP-violating correction to the dispersion relation at first order in gradients, and hence there is no 
CP- violating semiclassical force in the flow term of the kinetic equation at order H. This is in contrast 
to what was found in Refs. 0, H 0- 

To investigate whether there is CP- violation in the off-diagonal sector of the theory, we note that the 
mass matrix is diagonalized by the unitary matrix 



U 



cos 9 
sin 6 e la 



sin^e 
cos 6 



where 



tan 20 



2y\A*H 2 + f iH 1 
m 2 n — m 2 . 



tana 



with A*H 2 + fiHx 
easily obtain 



\A*H 2 + yU#i|e i<T , A* = \A\tf 



\A\ tan (3 sin a + \fi\ sin£ 
\A\ tan (3 cos a + \ fi\ cos( ' 

H = \fi\e^ and tan/3 = H 2 fH\. 



(4.38) 

(4.39) 
From this we 



iUd„lP 





— ie 



" T 



(1 - cos 26) sin2#e- i<T 
sin 26 e ia 1 - cos 26 



(4.40) 



Since the off-diagonal elements of S M are precisely the ones involved in the mixing of the diagonal and 
off-diagonal equations, H M is CP-violating whenever a ^ and tan/5 varies at the phase interface, 
which is in general realized at the phase transition in the MSSM by complex A and \x parameters. 
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5. KINETICS OF FERMIONS: TREE-LEVEL ANALYSIS 



In this section we consider the dynamics of fermions in the presence of scalar and pseudoscalar space- 
time dependent mass terms. Our analysis is of relevance, for example, for electroweak baryogenesis 
calculations, and establishes the nature of the propagating quasiparticle states in the electroweak 
plasma at the first order electroweak phase transition. In supersymmetric models baryogenesis is 
typically mediated by charginos and neutralinos j4^|. Since they both mix at the tree level, it 
is very important to consider the more general case of N mixing fermions J^. We derive the 
dispersion relation and kinetic equation for mixing fermions accurate to order h. We also construct the 
equilibrium Wigner function to order h, which can then be used for baryogenesis source calculations. 



This section is not just an overview of our original work |52j, |53|, but it also contains new insights 
and results. 



5.1. Spin conservation and the fermionic Wigner function 

In section l2~^l we have already derived the equation of motion (|2.55jl for the fermionic Wigner function. 
For the moment we neglect the interactions induced by C m t in the Lagrangean (J2HHJ). The equations 
for S K and S > are identical, so we omit the label in the following, indicating that the equations hold 
for both of them. The equation of motion is 

VS=(f+ % -$ -(m h + z 7 5 m a )e-5^ S = 0, (5.1) 

where we introduced the symbol T> for the kinetic derivative operator. For planar walls propagating 
in z-direction the mass terms in (|5.1j) simplify to mt l , a (z)e% 9zdkz , when written in the wall frame. 
Equation (|5.1jl is the tree-level fermionic master equation, which contains all the necessary informa- 
tion to study the nature of the propagating fermionic states in the presence of a space-time dependent 
mass term. We shall now show that, for planar walls, Eq. ()5.1|) contains a conserved quantity, corre- 
sponding to the spin pointing in the ^-direction in the rest frame of the particle. This will then allow 
us to recast S ,<,> in a form which is block diagonal in spin, such that different spin blocks completely 
decouple, provided, of course, they decouple at the boundaries. A very nontrivial consequence of this 
observation is that a quasiparticle picture, when formulated in terms of the spin states, survives to 
order h in gradient expansion. 

5.1.1. The 1+1 dimensional (k\\ = 0) frame 

In order to establish what precisely is the conserved quantity, we consider a particle with arbitrary 
momentum k and perform a boost into the frame where the particle momentum parallel to the wall 
vanishes, while the momentum component perpendicular to the wall remains unaffected: 

k = (k ,k x ,k y , k t ) -> k = (k ,0,0, k z ) , k = sign(k )Jk$ - k\ . (5.2) 
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In this frame the particle moves only perpendicular to the wall, which is why we call it "1+1 frame". 
The necessary boost is characterized by 



V|1 = F' 711 = r • (5 - 3) 

ko k 
In the spinor space it is represented by the operator 

L(k) = k ° + ~ k °-^-h (54) 
[2k (k + fc )]V2 

and its inverse L _1 (fc ,fc) = £(&0) — ^)> which act like 

L(k) f L-\k) = -i% - -f 3 k z = f . (5.5) 

The mass terms in the equation of motion are transformed as follows: 

L(k)m h , a (x)e-^ dk L-\k) = m h>a (x)exp( - % - d- (d~ k + [L{k)d k L~\k) , •]) ) . (5.6) 

Since L(k) is a function of ko and fcii, the commutator term does in general not even vanish for 
planar walls in the frames where m = m(t,z). The exception is the wall frame of a planar wall, 
where m = m(z), such that (dm) ■ (L(k)dkL~ l (k)) = 0, and the commutator in (J5.6|) vanishes. For 
this reason we work from now on in the wall frame of a planar wall. In section 15.2.41 below we 
remark on how to transform our results to other frames, in particular to the plasma frame. Upon 
transforming (|5.1jl . we then get 

VS = (i + - (m h + i^m a )^ dzdk ^j S = (1 + 1 dim. frame), (5.7) 

where the boosted Wigner function is 

S(k) = L{k)S{k)L-\k) . (5.8) 

Since the bubble wall, which is responsible for any dependence of the Wigner function on the average 
coordinate x^, is symmetric in the x-y plane, we can assume that the Wigner function has the same 
symmetry: 

S = S(k,t,z). (5.9) 

This form holds true for homogeneous boundary conditions, which is the case for most of the time 
during a first order phase transition, during which the bubbles are large, almost planar and far away 
from each other, when measured in the units of a typical diffusion scale. The derivative in ()5.7|) then 
reduces like fi — > 7°9f + 7 3 <9 2 , so that the spin operator in ^-direction 

S z = 7 °7 3 7 5 (5.10) 

commutes with the kinetic operator 

[V,S g ] = 0, (5.11) 
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and hence in this frame the spin in z-direction is a good quantum number. Assuming that the plasma 
is in thermal equilibrium before the phase transition takes place, the fermions are described by the 
equilibrium Wigner function 1)2.71]) . which is diagonal in spin in the 1+1 frame. The interaction with 
the bubble wall introduces no spin mixing, so that we can write S in the spin-block diagonal form 

S=J2$s, S S = P S SP S , (5.12) 



s=±l 

where P s denotes the spin projector: 

P s = ^(l + sS z ), P S P S , = 5 SS ,P S , s = ±l. (5.13) 

The spin diagonal Wigner function lives in the subalgebra of the the Clifford algebra of the Dirac 
matrices which is spanned by the matrices commuting with P s . We can write down a decomposition 
of S s by using a suitable hermitean basis of the subalgebra, 

S.(k) = iP.(k) [sj 3 j 5 g s (k) - sj 3 gl(k) + l^(fe) - tf%(k)] , (5.14) 

with g s a , a = 0, 1,2,3 being scalar functions. If we choose the following chiral representation of the 
Dirac matrices 

7" = I J",, °l ] , (5.15) 

where = (1, a 1 ), <x M = (1, — a 1 ), and a\ i = 1,2,3 are the usual 2x2 Pauli matrices, then the spin 
operator is 

S z = 7°7 3 7 5 = 1 ® ff 3 (5.16) 
and the spin block diagonality of the Wigner function becomes explicit: 

-il°S s = \p a ® (1 + 80*)? a . (5.17) 
5.1.2. The 3+1 dimensional (moving) frame 

Knowing that in the 1+1 frame the spin in ^-direction is conserved by the interaction with the bubble 
wall, we can construct the corresponding conserved quantity in the original frame by boosting the 
spin operator S z back to the original frame: 

S z (k) = L-\k)S z L(k) = i (k ol ° - fe„ ■ 7) 7 3 7 5 . (5.18) 

k v ' 

In order to clarify the physical meaning of this operator, we recall that the covariant form of the spin 
operator is given by the Pauli-Lubanski tensor 

S PL (k, n) = --p ^ 7 5 , e = (fe 2 ) 1/2 , (5.19) 
e 

which measures the spin of a particle with momentum k in the direction n (n 2 = — 1, n ■ k = 0). For 
simplicity we choose the normalization for Spl such that the eigenvalues are ±1. In literature one 
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often finds the normalization which corresponds to the eigenvalues ±h/2. In the on-shell limit, eo 

reduces to the particle's mass. In the rest frame the Pauli-Lubanski tensor becomes Spl — — > r y°rf'^f 5 , 
where n' is the spin vector in the rest frame, which is related to n by the corresponding Lorentz 
boost. To measure the spin in ^-direction in the 1+1-dimensional frame, we use 

& — (k , 0, 0, k z ), n" = -{k z ,0,0,k ), (5-20) 

e 

and we find, of course, SpL.(k,h) = 7°7 3 7 5 = S z . The same operator measures also the spin in 
z-direction in the rest frame, because the spin operator (|5.19|) is invariant under boosts as long as 
the direction of the spin vector is parallel to the momentum, n\\k. Now, by setting SpL(k,n) equal 
to S z (k) in (|5.18|) . we find that S z measures spin in the direction corresponding to 



n"(fc) 



k e 



( k k z \ 
k x k z 

k k 
ruyru z 

V ~ k o J 



(5.21) 



The same result is of course obtained by boosting n M in ()5.20j) to the original (moving) frame k^ . In 
the highly relativistic limit we have k\ —>■ k z , k$ — > k 2 , and our special spin vector (j5.21|) becomes 
proportional to k, such that the spin operator S z approaches the helicity operator, 

H(k) = -— ^7 5 
e 



k ■ 7°77 5 , h" = - U , (5.22) 




as one would expect. As usually, the helicity operator measures spin in the direction of a particle's 
motion, k = k/\k\. As a consequence, for light particles with momenta of order the temperature, 
k ~ T ^> m, the spin states we consider here can be approximated by the helicity states, which are 
often used in literature for baryogenesis calculations. 

The commutation relation ()5.11|) in the 1+1 frame now implies that by construction the spin operator 
S z (k) commutes with the kinetic differential operator T> in (|fl.ljl . provided the bubble wall is stationary 
and x-y-symmetric in the wall frame. We can slightly relax this condition: S z and T> commute even 
for non-stationary and x-y- dependent bubble walls, as long as the dependence is of the form 

S s = S s (k, t — v\\ ■ x\\, z) , (5.23) 

such that V|| = —vudt in D. So we can treat time dependent problems, which can be used to study 
how the system relaxes from some initial conditions to the stationary state, admittedly only for quite 
special forms of the initial conditions. 

As a consequence of the above discussion, even in the original moving frame, the problem splits into 
two non-mixing sectors labeled by the spin s: 

S=J2 S s, S s = P S SP S , (5.24) 

s=±l 



31 



where P s (k) = (t+sS z (k))/2, s — ±1 is the spin projector. We can again write down a decomposition 
of S s , 

S s = tP s [ S7 V<?o - *7 3 </S + IffJ - ^2] , (5.25) 

using a set of matrices that commute with P s {k). This is nothing else than the boosted 1+1 dimen- 
sional spin diagonal Wigner function: 

S s (k,x) = L(k)S s (k,x)L~ 1 (k) . (5.26) 

As a consequence of the hermiticity property ()2.51|) . the scalar functions g^(k,x) = gl(k,x) are all 
real. Since the direction in which spin is measured now depends on the momentum, it is not any 
more possible to find a representation of the Dirac algebra in which the block diagonality of S s is 
explicitly displayed. 



5.2. Constraint and kinetic equations 

We are now ready to study the tree-level fermion dynamics governed by Eq. ()5.1|) . In order to 
get the component equations, we take the block diagonal form (|5.25|) for the Wigner function S s , 
multiply ()5.1|) by P s (k){l, S7 3 7 5 , — is^ 3 , —7 s } and take the trace. The resulting equations are 

2ik g-Q - 2isk z g s z - 2im h g[ - 2irh a g s 2 = TrP s (k)lC^ (5.27) 

2ik Q g{-2sk z g s 2 - 2im h g s + 2m a g s 3 = TrP s (A;)s 7 3 7 5 ^ (5.28) 

2ikoti + 2sk z gt - 2m h gl - 2im a g s = TrP s (k) (-*s 7 3 ) (5.29) 

2ik g s 3 - 2%sk z gl + 2m h g s 2 - 2m a g{ = TiP s (k) (- 7 5 ) , (5.30) 

where we used the shorthand notations 

k = k + l - k ° dt — , k z = k z -~d z , (5.31) 

2 k 2 

and 

m hta = m h , a (z)e& 6 *> . (5.32) 

For completeness, we have added the contributions from the fermionic collision term. One should 
keep in mind that the dependence of the Wigner function on the time and parallel coordinates is 
restricted by (|5.23jl . 

In the general case with fermionic mixing, each of the functions g s a is a hermitean matrix in flavor 
space. The mixing is mediated through the off-diagonal elements of the mass terms ()5.32|) and thus 
appears already at the leading order in gradient expansion; hence it would be inappropriate to work 
in the interaction basis (j5.27M5.3THl without incorporating the flavor off-diagonal elements. Since when 
integrated over the momenta J d 4 k and f d A k{k^ /ko), equations (j5.27H5.3fHl yield fluid equations, the 
same conclusions hold for fluid equations. This is in discord with the strategy advocated in a recent 
work 51] on chargino-mediated baryogenesis in the Minimal Supersymmetric Standard Model, where 
it was argued that, when considering the dynamics of CP- violating sources, one should work in the 
weak interaction basis for charginos (without taking account of the off-diagonals). 
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Rather than considering the full dynamics of mixing fermions, we shall now argue that, in order to 
properly capture the dynamics of CP-violating densities to order h, it suffices to work in the spin 
and flavor diagonal bases, provided one transforms into the mass eigenbasis, in which m is diagonal. 
Here we generalize the analysis of Ref. (52] to 3+1 dimensions and include the space-time transients 
that accord with (|5.23p . More importantly, we fill in a gap in our original derivation in [52(1 . 



5.2.1. Flavor diagonalization 

In order to get to the mass eigenbasis, we diagonalize the fermionic mass matrix m. Since m is in 
general nonhermitean, the diagonalization is exacted by the biunitary transformation 

m d = UmV ] , (5.33) 

where U and V are the unitary matrices that diagonalize mm) and rwm, respectively. To make the 
analysis more transparent, we make an explicit separation of the spinor and flavor spaces by using a 
direct product notation ®. The master equation (|5.1|) then becomes 

(ft® 1 + S = C^, (5.34) 

where we have defined the mass term as 

m = Pr ® m + Pl ® rrv 

= 1 ® m h + ry 5 ® m a . (5.35) 
It is now a simple matter to see that the mass matrix m is diagonalized by the unitary matrices 

x = p L ® v + p r ® u = 1 ® hv + u) - 7 5 ® hy - U) 

Y = P L ® U + P R (8) V = 1 ® i(V + U) + 7 5 ® ~(V - U) (5.36) 

as follows: 

m d = XmY f . (5.37) 

The Wigner function then transforms as 

S d = YSX^ , (5.38) 

which is in general not diagonal. From here on we work in the frame where the bubble wall is at rest, 
so that we can make a spin-diagonal ansatz for the rotated Wigner function S s d which is analogous 
to (|5.25)1 . Note that the spin projector P s commutes with the rotation matrices. Since the mass 
term rh mixes spinor and flavor, the rotation matrices Y and X do so as well. As a consequence, the 
component functions g s ad of the rotated Wigner function are not just the flavor rotated components 
g s a of S s . Indeed, inserting ()5.25j) into ()5.38j) leads to the following relations 



9a* = \ [V(9o + 9l)V ] + U(g s - (5.39) 

9l d = \ [V(9s + 9l)V ] - U(g s - g°)U<<] (5.40) 

9[ d = \ [V{g{ - ig s 2 )tf + U(g{ + z^)yt] (5.41) 

9u = \ {Viti + i&tf + U (9 S 2 - ig{)V^] , (5.42) 
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that is <7g mixes with g| and g f mixes with g|. We now transform (J5.34j) by multiplying it from the 
left by X and from the right by to obtain 

® I + l -p-m d e& a **j S d = C H , (5.43) 

where C^ d = XC^Xt, and we defined a 'covariant' derivative 

D M = 1 ® Id^ - i[l ® E M , •] - ^{7 5 ® A M > "I > ( 5 - 44 ) 
where and A M are given in terms of the rotation matrices U and V as follows: 

A m = \ (V9^ - Ud^) (5.45) 

s M = ^ + ^c/t) . (5. 46 ) 

Note that the covariant derivative (|5.44jl can be obtained from 

<9 M m = <^(X+m d Y) = Xt(D^m d )Y (5.47) 

The higher order derivatives are obtained simply by iteration. Already from (|5. 43115. 4lj) it is clear 
that the kinetic and constraint equations, when written in the mass eigenbasis, will contain extra 
commutator and anticommutator terms. The choice of the rotation matrices U and V is not unique. 
After an x-dependent phase redefinition U — ► wU and V — > wV, where w is a diagonal matrix with 
eigenvalues of absolute value 1, U and V still diagonalize m. This freedom to redefine the rotation 
matrices was the source of some problems in finding the correct physical source in the WKB approach. 
We will find, however, that only the diagonal elements of the matrix A M are of relevance, and these 
are invariant under this reparametrization. 

We now project out the spinor structure in ()5.43|) by multiplying by P s (k){l, S7 3 7 5 , — is^ 3 , — 7 5 } and 
taking the spinorial traces. The procedure is identical to the derivation of (|5.27H5.3'0j) at the beginning 
of section 15.21 except for the subtlety related to the covariant derivative, because it contains spinor 
structure. The resulting equations are 

(2t~k - V-)g s M -s(2ik z + V z ^g s 3d - 2im hd e& 9 »>g' ld - 2tm ad e& 9 ^ g s 2d (5.48) 



TrlPC 



(2i~k - Vt)g{ d -s(2k z - iV^)g s 2d - 2im hd ei 5 ^ g s od + 2m ad e-^ A *g s 3d (5.49) 



Tr( S7 3 7 5 )P s C. 



ipd 



(2i~k - vf)g s 2d + s(2k z - iV z )g{ d - 2m M e& dk *g s 3d - 2im ad e^ 9 ^ g s od (5.50) 

= Tr (-is7 3 ) PgC^cL 

(2ik - Z>-)(& - s(2ik z + V^g'u + 2m hd e& 9 **g s 2d - 2m ad e^ 9 ^g s ld (5.51) 

= Tr(- 7 5 )P s C^, 

where we have defined the derivatives 

D z m hd = d z m hd - i[E z , m hd ] - {A z , m ad } 

D z m ad = d z m ad - i[T, z ,m ad ] + {A z ,m hd } (5.52) 
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and 



T>% = 7||3t + WW ■ V|| - zs[A„ •] (5.53) 
A + = 7||3t + 7p|| ■ V,| - is{A z , •} (5.54) 
V z = d z -i[Z z ,-} (5.55) 

Note that due to the 7 5 the anticommutator in ()5.44j) has become the commutator in ()5.53|) . Moreover, 
the derivative T>f is not hermitean. Indeed, the anticommutator term in (|5.54|) is antihermitean. 



5.2.2. Constraint equations 

The constraint equations correspond to the antihermitean parts of the equations ()5.48|l - ()5.5ip . Since 
we are interested in the dispersion relation accurate to first order in h, it suffices to consider these 
equations to first order in gradients: 

2k g s 0d - 2sk z g s 3d - {m hd , g s ld ] - - [D z m hd , d kz g s ld ] - {m ad , g s 2d } - - [D z m ad , d kz g s 2d ] = C s od 

(5.56) 

Zhgld + s{A z , g[ d } + s{d z - i[E z , -])g s 2d - {m hd , g s od ] - % - [D z m hd , d kz g s od ] (5.57) 

+^{D z m ad ,d kz g s 3d } -i[m ad ,g s 3d ] = C{ d 
2 ~ko9 s 2d + s{A z ,g s 2d } - s(d z - i[E z , -})g s ld - -{D z m hd ,d kz g s 3d } +i[m hd ,g s 3d ] (5.58) 

-{mad,g s od } - -[D z m ad ,d kz g s od ] = C s 2d 



Zkog s 3d -2sk z g s m + ]^{D z m hd , d kz g s 2d ) -i[m hd , g s 2d ]-^{D z m ad ,d kz g s ld } + i[m ad , g[ d ] = C s 3d , 



(5.59) 



where 



-4>d 



CS ° d = Y i Tll ( Ps(k)Clpd ~ P ° ik)C i 

C\ d = ^Tr(s^)^P s (k)C H -PKk)Cl d/ 
C* d = I T r(-^ 7 3 )(p s (A;)C^-Pt (A;)C t ( 

C s 3d = i Tr(- 7 5 ) (P.(*)C* - Pl{k)C\, d ) . (5.60) 

The trace is only to be taken in spinor space, not in the fermionic flavor space. For simplicity we 
now perform the analysis for two mixing fermions. Our findings are however valid for an arbitrary 
number of mixing fermions. Taking account of the fact that the off-diagonals are of the order H, 
equations (|5.56H5.59"j) now imply the following diagonal equations accurate to order H, 

2 ^o 9 s 0dn ~ 2m hdl g s ldn - 2m adl g s 2dn - 2sk z g s 3dn = C s odn (5.61) 

-2m M1 g s Qdn + 2~k g{ dn + 2sA zll g s ldu + sd z g s 2dn + {D z m ad ) u d kz g s 3dll = C[ dll (5.62) 

-2m adl g s odn - sd z g s ldll + 2~k g s 2du + 2sA zll g s 2du - (D z m hd ) u d kz g s 3dll = C s ldu (5.63) 

-2sA; < ^ 11 -(iP z m ad ) 11 9 fc ,^ dll + p z m M ) 11 9 fc ,y| dll + 2fco^Mii = C s 3dn , (5.64) 
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and similarly for the (22)-components. Solving the latter three equations in gradient expansion in 
terms of S^dii we find 



9idii ~ r 



1 

k 



#3<Z11 



1 

h 
l 

h 



m hdig s odll - isA zll (m wi ^ dll ) - - r sd z (m adl g s 0du ) 
k 2k 

s 1 
- rr(^ m ad) n (l + k z d kz )g s odu + -CJ^ 

m adi^odn + -^T sd z ( m hdi9 s mil ) - -~-sk zn (madig s M n) 
2k k 

s 1 
+ Tr{ D zm hd ) n {l + k z d kz )g s odu + -C s 2dll 
2k 2 



(5.65) 



(5.66) 



skz9odu + ■^- m hdi(D z m ad ) n d kz g s dn - -^-madiiD.nihd)^^^ + \c s 3dn ] .(5.67) 

2/Cn 2/Cn ^ J 



Here rrihdi denotes the first diagonal element of m^d, ""Wi is defined correspondingly. Inserting these 
relations into (j5.61|) we get the constraint for the diagonal densities 



tie - \m d \\ + f [\rm\KM* + 2A zii )}) g s odii 



k 

^Odii + 7~ { m hdiCldii + m adiC2dii + s ^zC 3dii ) 

ko 



(5.68) 



where we defined 



2 2,2 

,; = m hdi + m adi 



\m d \id z 6 di = m hdi d z m adi - m adi d z m hdi 



(5.69) 



where i = 1,2. Note that the energy shift in ()5.68j) can be written in terms of the rotation matrices 
U as 

\m d \l{d z e di + 2A 2il ) = -Z[U(md z mi)U% . (5.70) 



5.2.3. Quasiparticle picture and dispersion relation 

Equation ()5.68|) specifies the spectral properties of the plasma excitations. Since it is an algebraic 
constraint, the quasiparticle picture remains a valid description of the plasma to order h, as it was 
already pointed out in 0,0. Since the constraint equation contains no time dependence, it 
measures genuine spectral properties independent of the actual dynamical populations of the states. 
The homogeneous solution of the constraint equation (j5.68|) has the following spectral form 

g^ u (k,x) = 2n5(k 2 - Imatf + ^[Im^id^di + 2A zii )])\k \n sl {k,x) 
V k / 

= 2^ ^ 0=F ; ±8l) \~ko\n si (k,x), (5.71) 



with the following dispersion relations for fermions 



u si = u 0i - — [\m d \i (d z 6 di + 2A zii )] , (5.72) 

ZUJQiUJQi 
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where 



^oi = {k 2 + \m d \i)' 

wbi = {k 2 z + \m d ffi. (5.73) 



The normalization factors are 



7 s\m d \ 2 (d z 6 di + 2A zii ) 

Z si -l — . (5.74) 



J 0i 

We have normalized the solution ()5.71|) such that n si (k,x) represents the particle density in phase 
space {k, x}. We emphasize that the spectral solution (|5.71|) and the dispersion relation (|5.72|) are 
valid in general for all plasma excitations: equilibrium and stationary excitations, as well as for 
non-equilibrium, time dependent transients that conserve spin, that is for g$ d = g'Q d (k,t — k\\ ■ x\\, z), 
assuming of course planar symmetry in the wall frame. 

Since in the on-shell limit the Wigner functions S' < and S > satisfy identical Kadanoff-Baym equa- 
tions (|2.55p . and in addition they are related by the fermionic sum rule (12. 60(1 . based on the above 
analysis of S K we can easily reconstruct the spectral solution for the Wigner function S > , which we 
write as 

S> = £ ^s(k) - *7 3 ^ + ^9il ~ n 5 9%] , (5.75) 



=±i 

,>6 



where the spectral solution for g^ d reads 

= -2nJ2 6 [ k ° T y ±Sl) \ko\ [1 - n si (k,x)} . (5.76) 

Note that g^lu obeys the same equation ()5.68|) as gQ dii - Similarly, g^ d \ ± , g^ia an d 9s dii are related 
to gQ dii the same way as indicated in Eqs. ([5.65H5.67J) . 



5.2.4- Plasma frame 

An interesting question is of course how to generalize the solution for the Wigner function ()5.25j) to 
other Lorentz frames, an important example being the plasma rest frame. The natural solution is 
to apply a Lorentz boost on the Wigner function. The simplest boost operator corresponding to a 
boost v in z-direction can be easily constructed in analogy with (|5.4|) . and it reads 

w.) = 7+ P l 7 " + T)i^ - (577) 

where 7 = (1 - v 2 )- 1 ' 2 . The Wi gner function then transforms as 

S. = L-\A Z )S S L(A Z ) 

= -P s (k) [ 7 (1 + vyV) (a-yYSo - sj 3 g s s) + l?i ~ h 5 ~9i] , (5-78) 

where 

P s (k) = 7(fco r^ ) 7°7 3 7 5 " 7(1 + W)^ " 7||7 3 7 5 (5.79) 
k k 
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is the spin operator with k^ = 7 2 (/co — vk z ) 2 — fcjj, and g s a = g a (k,x). We leave as an exercise to the 
reader to find the direction of the spin vector corresponding to the spin operator (|5.79jl . To get the 
final form for the Wigner function in the new frame, one also needs to transform the coordinates and 
derivatives in (|5.65H5.67j) and (|5.71|) . which makes the final expression for iSf rather cumbersome. 
To get the Wigner function iS s in the plasma frame, one needs to choose the boost that corresponds 
to v = v w and 7 = j w = (1 — v"^)^ 1 ^ 2 . This exercise underlines clearly the advantages of workin g in 
the wall frame, in which the Wigner function has a particularly simple form. Some authors [50L l5ll| 
choose nevertheless to work in the plasma frame, although in a simplistic disguise. 

5.2.5. Kinetic equations 



We now take (minus) the hermitean part of Eqs. (|5. 48115. 5lf and get the following kinetic equat 



ions 



v t 9od + sV zg s 3d - {m hd sin (^D z d kz ), g[ d } + i[m hd cos (^D z d kz ), g{ d ] 

- {m ad sin (~D z d kz ) , g s 2d } + i[m ad cos {~D z d kz ),g s 2d \ = IC s od (5.80) 
( v t)h9 s i d + 2sk z g s 2d - {m hd sm (~D z d kz ), g s od } + i[m hd cos (-D z d kz ) , g s od ] 

- {m ad cos (-D z d kx ), gl d } -i [m^ sin {^D z d kz ),g^ d ] = K{ d (5.81) 
{ v i)h92d - ^sk z g[ d + {m hd cos {^D z d kz ), g s 3d } + i[m hd sin (~D Z d kz ), g s 3d ] 

- {m ad sin (~D z d kz ) , g% d } + i[m ad cos (~D z d kz ),g% d ] = JC 2d (5.82) 
v t9'ld + sV z9m ~ {m hd cos (-D z d kz ) , g s 2d } - i[m M sin (-D z d kz ),g 2d ] 

+ {m ad cos (-D z d k J, gl^ + iirriad sin (-D z d k z ), g[ d ] =)C s 3d , (5.83) 
where {T>f)h = J\\(dt — i[S t , ■]) + y\\v\\ ■ Vy denotes the hermitean part of the derivative, and 

KU = ~Trl(p s (k)C^ + PKk)Cl d ) 

*3d = -\^{si^){P s {k)C^ d + Pt{k)C\ d 

*Zd = —Tri-is^fah^ + P^kp^ 

Kid = ~\ Tr(-7 5 ) (Ps(k)C H + PKk)C\ d ) . (5.84) 

Again the trace must only be taken in spinor space. Since we are interested in the order h effects, we 
can truncate the kinetic equations (|5.8UH5.8'3*|) to second order in gradients. For example, the kinetic 
equation for particle density of spin s ()5.80|) reads 

(l\\d t + l\\v\\ ■ V|| - is[A z , -])g s 0d + s(d z - i[E z , -))g 
~ ^{D z m hd , d kz g{ d } + i[m hd , g{ d \ - - 



lid 



1 {D z m hd , d kz g{ d } + i[m hd , g{ d ] - \{D z m ad , d kz g 2d } + i[m ad , g 2d ] = K, s od , (5.85) 
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where we dropped the second order terms in the commutator, which is legitimate since they contribute 
only to the off-diagonal equations, and hence to order h 2 . 

From (J5.85)) we can immediately write the kinetic equation for the diagonal densities 
7|| (dt + #|| • V||)^ dll + sd z g s 3dn - (L> z m M )i 1 9 fc ^ 1 s dll - {D zmad ) n d kz g s 2dll 

T ^ s i^zl293d21 ~ ^*2103dl2 + ^zl29od21 ~ ^z2l9o d l 2 ) (5.86) 



1 



(D z m hd ) 12 d kz g s ld21 + (D z m hd ) 21 d kz g s ldl2 + (D z m ad ) 12 d kz g 2d21 + (D z m ad ) 21 d kz g 2dl2 



^Odll" 



In order to decouple the diagonal and off-diagonal equations we make use of the off-diagonal equations, 
which to leading order in gradients are 

i6(m hd )g s ldl2 + i5{m ad )g b 2dl2 (5.87) 
= -isA zl2 5{g s 0d ) + -{D z m hd ) 12 d kz Tr{g s ld ) + -(D z m ad ) 12 d kz Ti(g s 2d ) - isE zl2 5{g s 3d ) + K s Qdl2 

i6(m hd )g s 0dl2 + 2sk z g s 2dl2 - Tr(m ad )g s 3dl2 (5.88) 

1 i 
= -(D z m hd ) 12 d k Ti(g s 0d ) - -(D z m ad ) 12 d k J(g s 3d ) + K,\ dl2 

iS {m ad )g s 0dl2 - 2sk z g s ldl2 + Tr(m hd )g s 3dl2 (5.89) 

1 i 
= -{D z m ad ) l2 d kz Ti{g s Qd ) + -(D z m hd ) 12 d kz 5(g s 3d ) + lC s 2dl2 

Tr(m ad )g s ldl2 - Ti{m hd )g s 2dl2 (5.90) 

i 



% % 

-isE zl2 5(g s 0d ) + -{D z m ad ) 12 d k J{g[ d ) - -{D z m hd ) 12 d kz 5{g s 2d ) - isA zl2 S{g s 3d ) + K s 3dl2 



Combining ()5.87|) and ()5.90|) we obtain 

g s ui2 = flflj^ip j _ s ( D z m ad)i28(g s 0d ) - is[A zl2 5(m ad ) - iE zl2 Tr(m hd ))5(g s 3d ) 

l 

- -[Tr(m hd )(D z m hd ) 12 d kz Tr(gl d ) - Tr(m ad )(D z m ad ) l2 d kz 5(g s ld )\ 
~ ^[^(m hd )(D z m ad )i 2 d kz Ti(g 2d ) + Tr(m ad )(D z m hd ) 12 d kz 5(g 2d )} 

- iTr(m hd )JC s 0dl2 + 8{m ad )lC s 3dx ^ (5.91) 
9 S 2di2 = s(\m d \) 2 { s( ^ Dzmhd ^ 125 ^d) + is[A zl2 8{m hd ) + iE zl2 Tr(m ad )]8(g^ d ) 

- ~[Tr(m ad )(D z m hd ) 12 d kz Tr(gl d ) + Tr(m hd )(D z m ad ) 12 d k J(gl d )} 

- -lTr(m ad )(D z m ad ) 12 d kz Tr(g 2d ) - Tr(m hd )(D z m hd ) 12 d kz 8(g 2d )} 

- iTr( mod ))C s 0dl2 - 5{m hd )K s 3dl2 ) , (5.92) 



where 



S(\m d \) 2 = Tr(m hd )5(m hd ) + Ti(m ad )S(m ad ) . (5.93) 
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From (I5.871l5.90j) and (151)11151)21) we then get 



S&fta = j^^{-2k z [i: zl2 5(g s od ) + A zl2 5(g^ 



+ sk z [(D z m ad ) l2 d kz 5(g s ld ) - (D z m hd ) 12 d k J{g 2d )} 

- -[Tr(m hd )(D z m a d)i2 - TT(m ad )(D z m hd )i 2 ]d kz S(g 3d ) 



- iTr(m hd ))C s ldl2 - iTi(m ad )}C s 2dl2 - 2isk z JC s 3dl2 j (5.94) 

9 s 3dl2 = jr^{-2k z [A zl2 5(g s J + E zl2 6^ 

- isk z {(D z m hd ) 12 d kz Tr(g s ld ) + (D z m ad ) 12 d kz Ti(g 2d )} 
+ -[5{m ad )(D z m ad ) 12 + 5{m hd )(D z m hd ) 12 ]d kz 5{g s 3d ) 

- 2isk z K s QdX2 - 5{m ad )K s ldl2 - 5(m hd )lC s 2dl2 j . (5.95) 

Equations (j5.91H5.95j) are the off-diagonal densities correct to leading (first) order in gradients. We are 
now ready to compute the contribution from the off-diagonal densities in the kinetic equation (j5.86j) . 
Makine use of the hermiticity properties A Z J 2 — A z2 i, — ^z21j 9adi2 — 9ad2ii 

and (15.91115.951) . 

after some algebra one finds that the contribution from the off-diagonal densities vanishes: 

^'zl25 , 3d21 ~~ ^z21#3dl2 + ^zl29od21 ~ ^z2\9odV2 = 

(D z m hd ) 12 d kz g[ d21 + (D z m hd ) 21 d kz g s ldl2 + (D z m ad ) l2 d kz g 2d21 + (D z m ad ) 21 d kz g 2dl2 = , (5.96) 

where we neglected the collision terms that are suppressed by derivatives. With the help of this 
remarkable result and the constraint equations (j5.65H5.67j) the kinetic equation (j5.86j) finally reduces 
to the familiar form 0, E3j 
1 

k 

Note that this equation holds for both g^ d and g^ d . Recall furthermore that for spin conservation 
the space-time transients are constrained to satisfy gQ dii = go dii (k,t — v\\ ■ x\\,z). The collisional 
contributions from the constraint equations are higher order in gradients and can be consistently 
neglected. 

The functions g s ad (a = 0, 1, 2, 3) satisfy four apparently different equations (j5.8Uf!5~5Hj) . These func- 
tions are related by the constraint equations (j5.56ti5"59j) . which reduce the number of independent 
functions to a single one, g$ d , projected on the quasiparticle shell (j5.68j) . We will show in an appendix 
in Paper II that all four kinetic equations are actually mutually dependent in a self-consistent man- 
ner, and thus equivalent to the kinetic equation for g^ d , plus the on-shell equation (j5.68j) . Our proof 
includes not only the flow term, but also the collisional sources. The equivalence of these equations 
including the collision term is a very nontrivial consistency check of our approach to the kinetics of 
fermions. 

5.3. Boltzmann transport equation for CP-violating fermionic densities 

Since our primary interest is transport of CP- violating densities, we show here how - starting with the 
Kadanoff-Baym equations for fermions (j2.55j) (cf. also (j5.97j0 - one obtains the on-shell Boltzmann 



k-d- ^(d z \m d \* - td z [\m d \Ud z 9 dl + 2A zli )] )d kz ]g s odii = lC* dii . (5.97) 
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transport equations. The crucial difference with respect to the scalar case discussed in section 14. f I is 
the CP- violating semiclassical force in the flow term of the kinetic equation for fermions. The force 
appears at second order in gradients, or equivalently at first order in an expansion in h. 

We begin our analysis with recalling the C and CP transformations for the fermionic fields, 

i) c {u) = Cip(u)tf = i7°7 2 V> T (w), 

V> c (w) = Ci){u)tf = ^ t (m)z 7 V , (5.98) 

and for Dirac's 7-matrices, 

cy T ct = * 7 ° 7 v T (* 7 V)t = - 7 M 

C 7 5 C f = i7 7 2 7 5 (i 7 7 2 ) t = 7 5 . (5.99) 

In the relativistic limit the Standard Model fermions are chiral and couple differently to the weak 
interactions, so that both charge C and parity P symmetries are strongly violated. But the combined 
symmetry CP is violated only very weakly, hence it is natural to consider antiparticles to be related 
to particles by a CP transformation. Therefore we also quote the parity transformations, 

^{u) = Vip(u)V^ = 7 ( V(«), TVPt = 7 V 7 ° = y 4 = 7m 

tfj p (u) = Vi){u)V^ = tp(uh° , P 7 5 7 3t = 7 ° 7 5 7° = -7 5 • (5.100) 

Combining these relations with the definitions of the fermionic Wightman functions (|2.5j) . we easily 
find the relations for their C- and CP-transformations: 

S K (u,v) 7 ° 7 2 ^ >T (^,«)7V (5.101) 
S < (u,v) ^ - 7 2 S >T (t;,M) 7 2 , (5.102) 

which in the Wigner representation become 

S < (k,x) 7 ° 7 2 ,S >T (-fc,x) 7 7 2 = S c< (k,x) (5.103) 
S < (k,x) ^> - 7 2 5 >T (-fc,x) 7 2 = S cp< (k,x); (5.104) 

analogous relations hold for S c> and S cp> . As we did in the scalar case in (j4.24ff4~23)) . we employ 
an additional inversion of the spatial parts of position x M = (x°, — x l ) and momentum k^ = (k°, —k l ) 
in the definition of S cp . In the case of several mixing fermions, the Wigner function is a matrix in 
flavor space which is affected by the transposition, too. 

We have argued in section |3] that a weak coupling reduction to the on-shell limit of the equation of 
motion ()2.55|) for the fermionic Wigner function results in the equation 



ft+-0 x - (m h (x) + i^m a (x))e~^ dk ^S > ' < (k,x) =C i ,(k,x). (5.105) 
Taking account of the hermiticity property ()2.55|) . hermitean conjugation and transposition lead to 
f ~\€ ~ KW +tl 5 m* a {x))e^* dk ^S > ' < {k,x) T = - 7 °C;(A;,x) 7 . (5.106) 
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By commuting 7°7 2 through and reversing the sign of the 4-momentum, we get the following equation 
for the C-conjugate of the Wigner function ()5.10Hj) . 

( f + l - I - (ml{x) + i^m^x)) e-^ 9 *)s c <>>(k, x) = i 2 C;(-k, x)j 2 . (5.107) 

Similarly, by commuting r y 2 through ()5.106|) . we arrive at the equation for the CP-conjugate Wigner 
function (I5.104|) : 

K(*0 - rf<( x )) ^ %9 ^S c ^>(k,x) = -7°7 2 C;(-M)7°7 2 • (5.108) 

Several remarks are now in order. Note first that the presence of imaginary elements in the mass 
matrices, rrih = (1/2) (m + m)) and m a = (l/2i)(m — m^), may violate both C and CP symmetry. 
This is true, of course, provided the imaginary parts cannot be removed by field redefinitions. Second, 
while C is not in general violated by a (real) antihermitean mass term m a , CP is violated, provided 
m a = m a (x) is space-time dependent, such that it cannot be removed by field redefinitions jsil ]. 
Finally, Eqs. ()5.107H5.108|) illustrate where in the collision term are potential sources of C and CP- 
violation. As expected, both C and CP may be violated in the presence of complex Yukawa couplings, 
appearing in collision or mass terms. This can be seen from the complex conjugation of the collision 
term, C^. Further, the weak interaction vertices contain Pl,r, which is also a property of our model 
Lagrangean ()2.14j) . Now since 7 2 -Pl,_r7 2 = Pr,l flips chirality, while (7°7 2 ).PL,i?(7 7 2 ) = Pl,r leaves 
it invariant, we conclude that the weak interaction vertices violate C, but they are invariant under 
CP, as they should. Using the explicit expressions for the collision terms in Paper II and ()5.1.2|) . 
one can show that the right hand side of equation ()5.108|) is equal to C^, with all Wigner functions, 
both scalar and fermionic, replaced by their CP-conjugate counterparts, and the Yukawa coupling 
matrices y replaced by y* . Further details of the analysis of the collision term are left for Paper II. 

A careful look at equation ()5.108|) shows that the kinetic operator for S cp {k) commutes with P s (k) 
rather than with P s {k). So the spin-diagonal ansatz for the CP-conjugate Wigner function, already 
written in the mass diagonal basis, is 

s s 

Now there are several ways to obtain the constraint and kinetic equation for the component function 
fi'odii- First, one could repeat all the steps from the previous section. Since the mass term in the CP- 
conjugate equation (|5.108j) is the complex conjugate of the mass term appearing in the equation for 
the original Wigner function ()5.105j) . we would begin by carrying out the flavor transformation, this 
time using the rotation matrices X* and Y*, then taking the appropriate spinorial traces. The anti- 
hermitean and hermitean parts would lead to constraint and kinetic equations for the CP-conjugate 
component functions, respectively, and finally to the desired equations for g^ii- 
Another possibility is to note that equation ()5.108|) for S cp is obtained from equation (|5.105|) for S by 
replacing the mass by its complex conjugate and replacing all explicit occurrences of the momentum 
k by k, as well as sending all spatial derivatives d x to d s , and finally replacing the collision term on 
the right hand side, as indicated in Eq. ()5.108j) . When we apply exactly the same changes to the 
constraint and the kinetic equation for g d, which includes 9^—9 and A z — > —A z because of the 
complex conjugated mass, we find the respective equations for g c ^ d . 
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The simplest method, however, is based on the observation that, because of relation (J5.1U4j) . the 
Wigner function S contains complete information about its CP-conjugate S cp , and the same then 
holds for the component functions go and g^ as well. Relation 1)5.104)1 is covariant under the flavor 
rotation, 

Sj p< (k,x) = Y*S cp< (k,x)X. T = - 7 2 (YS > (-k 1 x)X^) T l 2 = ~l 2 S> T (-k, x) 7 2 , (5.110) 
so we have to match the decomposition (J5.109)) for S cp< with 

- 7 2 [£>,(-£;, x)) V = £ iP.(-k) [ ~ *tV<&[? + s 7 3 ^7- + lg>£ + vfg>£] (-*, x) . (5.111) 

s 

This implies the identification 

9%Z(k,x) = -g>£(-k,x) i (5.112) 

where we used P s (—k) = P B (k). Looking at the constraint equation ()5.68)) or the spectral ansatz ()5.76|) 
for <7ota i; we immediately see that the component functions of the CP-conjugate Wigner function live 
on the same energy shells u si as the ones of the original Wigner function. In other words, the 
dispersion relation for antifermions is identical to that for fermions, a conclusion that can also be 
reached by considering the CP-conjugated constraint equations. Therefore we can make the spectral 
ansatz 

g<£(k, x) = 2nJ2 6 [ k ° T ^ ] \k \ n%(k, x) . (5.113) 

Inserting this and the spectral ansatz (|5.76|) for into ()5.112|) shows that the particle density at 
negative energies is related to the density of antiparticles. For k on the positive or negative energy 
shell we have 

n^(k , k,x) = 1 - n^ si ( - k , -k,x) . (5.114) 

In order to find the kinetic equation for the CP-conjugate component function, let us first recall the 
kinetic equation ()5.97)) for gQ dii from section IB. 2. 51 

' \k ■ d- \{d z \m d \l - td z [\m d \*(d z 6 dl + 2A zii )] )d k ]g s odii (k, x) = JC s odii (k,x) . (5.115) 
- v k 



k 



Again we make use of relation ()5.112j) and find 

l\k ■ d- Ud z \m d \* - f d z l\m d \1(d z 6 dl + 2A zii )] )dJgZl i (k,x) = -K^±(-k,x) . (5.116) 
ko L i \ k > J 

Needless to say, all ways of obtaining these equations yield the same results, although this might not 
be obvious for the collision term. Starting with the right hand side of equation (|5.108|) . one would 
expect to find 

/C^(A;,x)EEiTr(lP s (%Vq dii (-A;,x) 7 V + h.c.) (5.117) 

as the collisional contribution to the kinetic equation for g c Q v . But since this is a trace in spinor space 
and we are looking at the diagonal elements in flavor space, we can apply a transposition to find 



= ^Tr(( 7 VP s (%V) T C; dii (-A ; ,x)+h.c: 
= ^Tr (p_ s (fc)C^ ii (-fc,x)+Pl s (A;)^4 i (-fc,x)) = -IC ^(-k,x) . (5.118) 
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We define the distribution function for particles as the projection of the phase space densities on the 
positive mass shell, 

fai+{k, x) = n si (u si , k, x) , (5.119) 
and then the kinetic equation for c/q P suggests to define the distribution function for antiparticles as 

f si -{k,x) = ri%(u} s i, k,x) = 1 - n^ si (-uj si , —k,x) . (5.120) 

We will verify in Paper II that f S i±(k) indeed measure the density of particles and antiparticles 
with momentum k and spin s. To obtain the Boltzmann equations for these densities we have to 
integrate ()5.115|) and (|5.116|) over positive frequencies: 

1 f°° dk 

^(9 i + i'V J + FA)/«+= / — ^ dii (M (5.121) 

^si Jo 71 

1 f°° dk 

-(^i-Vi + FA)/,^ / —ICfiufax), (5.122) 

^si Jo 71 

with the quasiparticle velocity and semiclassical force 



k 



(5.123) 

F si = -J-fo| m< ,|? - ^d z [\m d \*(d z 9 di + 2A 2ii )] ) . (5.124) 

Note that, in addition to the classical terms, both v S i and F S i contain spin dependent contributions 
arising from derivatives of the pseudoscalar mass term, and thus a potential for inducing CP- violating 
effects. 

The Boltzmann equation ()5.122|) for the antiparticle densities can also be obtained by integrating the 
kinetic equation for g over negative energies, and then using (|5.120|1 . This procedure has often been 
advocated as a way of identifying CP-violating densities. It works, because the complete information 
about S <cp is contained in S > , which in the on-shell limit is completely determined by S K , as can be 
shown with the help of the fermionic spectral sum rule (|2.60|) . 

For practical baryogenesis calculations one often deals with problems close to thermal equilibrium, 
in which case it is convenient to treat the Boltzmann equations (J5.121j) and ()5.122|) in linear response 
approximation with respect to deviations from thermal equilibrium. We note that the thermal 
equilibrium distribution function in the wall frame has the form 

^ = e P-y w (u si +v w k z ) _|_ (5.125) 

with (3 = 1/T, v w = v w z is the wall velocity and ^ w = (1 — ^) -1 ^ 2 is the corresponding boost factor. 
This is a generalization of the thermal equilibrium distribution function (|2.73J) which includes both 
the effects of CP- violation from a pseudoscalar mass and of a moving plasma on thermal equilibrium. 
We show in Paper II that this form for the equilibrium distribution function leads to the correct 
energy conservation in the collision term. Furthermore, the distribution function ()5.125j) satisfies 

Z^(d t + v si ■ V, + F si d hz )f:? = -v w f£(l - ft)2tfPF« (5.126) 
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such that, for a wall at rest, Eq. (|5.125|) satisfies the flow equation exactly, as one would expect 
from the correct equilibrium distribution function. Note that the degeneracy between the states with 
opposite spin is broken by this equilibrium, because particles with opposite spin satisfy different 
energy-momentum dispersion relations. 

Our main interest is the transport equation for a CP-violating distribution function that can even- 
tually lead to the creation of a Baryon asymmetry. We first define 

fsi± = m + l5fr n + 5f si± , (5.127) 

where 5f S i± are spin dependent correction functions. Since spin dependence occurs only at order h 
in our Boltzmann equations and all sources are wall velocity suppressed, as for example can be seen 
in (j5.126|) . we find 5f si ± = 0(v w h). The distributions f si ± can also contain a CP-even deviation from 



thermal equilibrium, $f? ven , which is relevant for the dynamics of phase interfaces [40, 74j,|75j. Note 
that, since Sff ven is not suppressed by powers of h, it does not in general decouple from the equation 
for the densities Sf S i±. In order to properly treat the dynamics of CP-violating densities, one has to 
solve also for the dynamics of CP-conserving densities! Together with Sf^ ven = 0(v w h°), Eq. (j5.126J) 
implies that it suffices to solve the equation for Sf^ ven to leading (classical) order in gradients, which 
is obtained by adding (I5.121|) and (I5.122jl . 

F) \m I 2 f°° rlk 

(d t + v 0i ■ V s + F 0t d kz )Sfr n + (3 foil - /oK^ 1 ^ = / — (JC s 0di± (k) + K^{k)) , (5.128) 
where 

voi = — , ^oi = 7; (5.129) 

u) i ZujQi 

are the classical particle velocity and the classical force, respectively, and f 0i = l/(e^ 0i + 1). 
Based on the first order correction 5f S i±, we can form two CP- violating distribution functions: 

5f si = 6f sl+ -6f sl . (5.130) 

5f a si EE rt/,;. ~df ,, . (5.131) 

As we will explain in Paper II, Sf^ and 5f^ are related to the vector and axial vector density in 
phase space, respectively. Working to second order in gradients and to linear order in v w , we obtain 
the equation for Sf^ by subtracting (j5.122j) from (|5.121j) : 

{d t + v oi ■ V, + FudkJSfc = / ^ (^odii(^) - ^aiW) • (5-132) 

Jo 71 

Note that the flow term contains no source whatsoever. This means that a CP-violating semiclassical 
force (J5.124j) cannot create vector charge densities (this remains no longer true, however, when one 
includes the possibility of coherent particle production, which can be exacted by a dynamical tracing 
of flavor mixing [7^). 

The transport equation for df^ is obtained by changing the sign of s in ()5.122|) and then subtracting 
it from (15.1211) . With the help of (I5.12(ij) we find 

r°° rlk 

{d t + v 0i ■ V, + F 0t d kz )6f: t + sS?™ = / {}C s odii (k) - fCZim) . (5.133) 

Jo 71 
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This equation now has a source: 

Sf ow = -lv w f 0l (l-f 0l )(3\6F i + F 0l 6 -^ + F 0t 5u Jl (3(l-2f 0t )+F 0t 6Z l 



+ 



-\5Zi(d t + v 0i ■ V a? + F 0i d kz ) H c n , 



F ^ + ^ 



5 fc2 \Sf t 



(5.134) 
(5.135) 



where the definitions 



8uji 



d z [\m d \ 2 (d z 9 dt + 2A zii )} 
[\m d \f{d z e di + 2A zii )] 



5Zi = SiOi 



(5.136) 



have been used. We now pause to comment on the physical meaning of the various terms appearing 
in the Boltzmann transport equation (|5.133|) . First, the CP- violating density, Sf^, evolves on phase 
space according to the standard Boltzmann flow derivative, d/dt = d t + v 0i ■ + F 0i dk z - Second, 
various CP- violating sources are displayed in ([5.134)1 . The first term, 8Fi, is the source arising from 
the semiclassical force that has been usually taken into account in the WKB-literature 

0Q- Note 

that the second term has a similar origin. Indeed, it arises from the CP-violating deviation in the 
energy-momentum relation appearing in the force term, — (d z \m d \i)/2u S i. The third term originates 
from the CP- violating split in the dispersion relation in the equilibrium solution (|5. 125(1 . and it thus 
vaguely resembles a local version of spontaneous baryogenesis. And finally, the fourth term arises 
from the gradient "renormalization" of the Wigner function. All four source terms are of the same 
order in gradients, and hence, a priori, they are all equally important. In Paper II we take moments 
of the Boltzmann equation 1)5.133(1 in order to obtain fluid equations. Since the source term 1)5.134)1 
is symmetric in momentum, it will contribute to the zeroth moment equation. The parametric 
dependence of this source becomes explicit by rewriting it like 



d 3 k 



(2n] 



cflow 

3** 



d z [\m d \*(d z 9 di + 2A zii )} 



(2 



7T 



+ f3 2 {d z \m d \ 2 t ) [\m d \ 2 (d z 6 dl + 2A zii )] 



I,, 



(5.137) 



In figure El we plot the dimensionless integrals X a and as a function of the mass. When \m d \i T 
(|^d|i 3> T), the source X a (2j) dominates. The sources are comparable in strength when \m d \i ~ 
T. Third, the terms ((5.135(1 that couple the CP-even deviation form equilibrium have never been 
considered in literature. Nevertheless, they are formally of the same order as the source terms ((5.134(1 . 
and thus cannot be neglected. While the first term in ((5.135(1 can be reexpressed in terms of the sum 
of the collision terms, as indicated by Eq. ((5.128(1 . the latter terms cannot, and represent genuine 
dynamical source terms, which have so far not been included in baryogenesis calculations. Fourth, 
the right-hand-side of Eq. ((5.133(1 contains the collision terms, to be considered in detail in Paper II. 
Finally, we remind the reader that all of the source terms in Eq. ((5.133(1 are second order in derivatives 
(first order in H), and are suppressed by the wall velocity. 



5.4. Applications 

In this section we have developed a formalism for the treatment of CP-violating effects of space- 
time dependent scalar and pseudoscalar mass terms in kinetic theory. We have shown how the 



46 




(3\m 



FIG. 5: The integrals Z a and Z& in Eq. (|5.137j) as a function of the mass. We scaled I a with \m\ 2 because this 
is the way it appears in the source term, analogously If, is scaled with |m| 4 . Note that these contributions 
enter the flow source with different signs. 

effect of a CP-violating shift in the dispersion relation induces an order h semiclassical force in the 
flow term of the Boltzmann transport equation. The force acts on particles and antiparticles of 
opposite spin in opposite directions. Since we have included the possibility of fermionic mixing, this 
formalism is suitable for studies of baryogenesis in supersymmetric theories at a strongly first order 
electroweak phase transition, in which baryon production is typically biased by CP-violating chargino 
or neutralino currents. Further, our formalism is also suitable for studying baryogenesis problems 
from coupling of fermions to the bubble wall in two Higgs doublet models. For completeness we shall 
now discuss in some detail the case of charginos in the Minimal Supersymmetric Standard Model 
(MSSM) and its nonminimal extension, and finally remark on baryogenesis sources in two Higgs 
doublet models. 



5.4-1- Minimal Supersymmetric Standard Model (MSSM) 

Now that we have a general expression for the semiclassical source in the case of mixing fermions, 
we want to study two explicit examples, which are of relevance for baryogenesis. First we compute 
the source in the transport equations for the chargino sector of the MSSM. The chargino mass term 
reads 

ip R mip L + h.c. , (5.138) 
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where ipR = {W^, hf R ) T and ipL = (Wti ^2t) T are the chiral fields in the basis of winos. The mass 
matrix is 

where H\ and i?2 are the Higgs field vacuum expectation values and /x and m 2 are the soft super- 
symmetry breaking parameters, which introduce CP- violation 0]. Since for a reasonable choice of 
parameters there is no transitional CP-violation in the MSSM, we can take the Higgs expectation 
values Hi and H 2 to be real j40, 77]. The matrix that diagonalizes mm^ can be parametrized as ji^] 



U= - ^ f^< A+ , A > ° V (5.140) 

where 

a = gijn^Hi + fi*H 2 ) 

A = | m2 |2_| / ,|2 + ^2_^2 ) 

A = V^ 2 + 4|a| 2 , (5.141) 
and /ij = \Hi\. The mass eigenvalues of the charginos are given by 

ml = \ (M 2 + |/f + g\h\ + hi)) ± I . (5.142) 

Upon inserting (j5.139M5.142j) into (j5.70j) . it is straightforward to show that the chargino source term 
can be recast as 2 

[\m d \ 2 (d z e + 2A z )] ± = Tj^m 2 )d z (hih 2 ) . (5.143) 

The sources figuring in the transport equation written for charginos in (j5.134j) can be easily recon- 
structed from equations (I5.141j) and (j5.143j) . The result (j5.143j) agrees with the one found by WKB 
methods in j^J. In however a different dependence on the Higgs fields was obtained. 

5.4-2. Nonminimal Super ■symmetric Standard Model (NMSSM) 

We now consider an extension of the MSSM which contains a singlet field S in the Higgs sector, which 
can induce additional CP-violation in the Higgs sector. In particular, the Higgs vacuum expectation 
values may become complex. The singlet couples to higgsinos, and therefore we obtain the mass 
matrix by generalizing the higgsino-higgsino component of the chargino mass matrix (j5.139j) 

fi — > jx = /i + XS , (5.144) 

where A is the coupling for the higgsino-higgsino-singlet interaction. The field content we consider is 
the same as in the MSSM, so the mass matrix is 



m= I m2 9m I . (5.145) 

1 gH{ p, 1 
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This matrix is still diagonalized by U in (J5.14U)) . We write the Higgs expectation values as 

H i = h i e idl , i = l,2, (5.146) 
where only one phase is physical. With the gauge constraint [3| 

h\Q\ = h\6' 2 (5.147) 

we can write 

e[ = f 2 e' , 0' 2 = §0', (5.148) 

where 6 = 61 + 62 is the physical CP- violating phase, and h 2 = h 2 +h?>. Now everything is prepared to 
write the NMSSM-source term. It can be divided into three contributions, which have to be added. 
The first one is a generalization of the chargino source (|5.143jl 



\m d \ 



(d z 6 d + 2A z )) hih2± = T^ifim^ihh) 1 (5.149) 



for the case involving a new scalar field S and complex Higgs expectation values. In addition to this 
there are two new types of sources. One of them is proportional to a derivative of the CP-violating 
phase 6 in the Higgs sector: 

(\m d \ 2 (d z 6 d + 2A z )) e± = ((A ± (|m 2 | 2 + T ^m 2 e i6 )h x h^ . (5.150) 



Finally, there is a source that can be written as a derivative of the singlet condensate: 

) =±^ 



{\m d \ 2 (d z 6 d + 2A z )) s± = ±±f$t(m 2 H 1 H 2 S') (5.151) 



Xg 2 
H — — 
2A 



(A±(|/2| 2 + ^ 2 -|m 2 | 2 ))s(/i^'). 



The mass eigenvalues m±, that is the diagonal elements of |m,i| 2 , can be obtained from the corre- 
sponding expression f)5. 142j) in the MSSM part with the replacement /j, — > fl. 
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6. CONCLUSION 



In this work we perform a controlled first principle derivation of transport equations for a model 
Lagrangean of chiral fermions Yukawa-coupled to a complex scalar. Our treatment is accurate to 
first order in an h expansion, or, more concretely, to first order in a gradient expansion with respect to 
a slowly varying scalar background field, which is formally valid when hd <C Hk. Being valid to order 
h, our treatment allows us to trace the propagation of CP-violating fluxes, which is of relevance, for 
example, for electroweak baryogenesis at a first order phase transition. We consistently include the 
collision term, although in this first paper of our series only in a formal way. The actual evaluation 
of these terms can be found in Paper II. 

We address an open question of electroweak scale baryogenesis mediated by mixing fermions, which 
couple in a CP-violating manner to a propagating bubble wall of a first order phase transition: which 
basis should be used to model the kinetics of mixing fermions? We show that, if one is limited to a 
diagonal approximation, the mass eigenbasis is singled out as the only basis in which the diagonal and 
off-diagonal elements of the distribution function decouple at order h in a derivative expansion. Our 
derivation is valid when there are no nearly degenerate mass eigenvalues, that is when hk-d <C 5(m^), 
where S(m^) denotes the (minimum) split in the mass eigenvalues. No such claims can be made for the 
flavor (weak interaction) basis, in which flavor mixing is present already at the classical level 0(h°). 
This indicates that the use of the flavor basis in transport equations is at best problematic, unless 
flavor mixing is consistently included. Of course, the final resolution of this problem can only come 
from a basis independent treatment, which would include the dynamics of both flavor off-diagonal 
and diagonal CP-violating densities. At the moment no such treatment is available, however. 

We also address various other issues, which comprise a proof that the Kadanoff-Baym equation for 
a single scalar field can be reduced, in a weakly coupled regime, to an on-shell Boltzmann equation, 
which includes both the self-energy and the collision term, approximated at the same order in a 
coupling constant expansion. An inclusion of the hermitean self-energies would be desirable, but 
has not been done yet. These self-energies mix spin and moreover provide an additional source of 
CP-violation, which we expect to be of a similar strength as the source in the collision term (see 
Paper II). Further, we demonstrate that no CP- violating source is present at order h in the flow 
term of the scalar kinetic equation. We then derive a Boltzmann transport equation for the relevant 
fermionic quasiparticles with a definite spin. We include all of the CP-violating sources arising from 
the flow term present at the order h, which, of course, include the semiclassical force, but also a few, 
as-of-yet unencountered, sources. One of the new sources is induced by a CP-even deviation from 
the equilibrium distribution function, which itself is formally of order h°. 

Our method represents a formalized and controlled implementation of the original heuristic semi- 
classical (WKB) treatment of the problem 34,|35j], in which one calculates the relevant (fermionic) 
dipersion relation accurate to order h by the means of a WKB analysis, and then inserts the result 
in the appropriate kinetic equation. As we have shown here, when applied with due care, the semi- 
classical approach leads to a correct semiclassical force in the flow term, but it does not permit the 
treatment of collisional sources, nor does it allow a critical assessment of the quasiparticle picture 
of the plasma or the implementation of any effects beyond the quasiparticle picture. None of these 
limitations apply for the approach presented here. 
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Our work is completed in Paper II. There we discuss the collision term, which in the present paper 
was formally maintained in all equations, but never evaluated explicitly, and we derive and study a 
set of fluid equations, which form an approximation to the Boltzmann equations found here. 
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APPENDIX A: GRADIENT EXPANSION IN THE OFF-DIAGONAL SCALAR KINETIC 
EQUATION 



We shall now prove that, provided k ■ d, (<9M|) • d k <C <5(Mj), gradient expansion applies for the 
off-diagonal equation. Equation (|4.1 6)1 is of the form 



k ■ d + l -{8Ml) ■ d k + l -5{Ml) - ik ■ 5(E) )zA< 



S 



12- 



(Al) 



where Si 2 = 5i 2 [A d ] represents the source composed of the diagonal elements. This equation can 
be solved by using the method of Green functions as follows. Consider for simplicity the following 
problem 

(k-d + ~5(M 2 ) + e)G r e (k;x,x') = 5\x-x r ), (A2) 
which is solved by the retarded Green function 



G[{k- x, x') = ^-6(t - t')5 3 (f - f - if) ) exp ( - 

k \ k 



2 k n 



+ e 



(*-0): 



(A3) 



where e represents an (infinitesimal) positive dissipation term. Formally, this solution contains rapid 
oscillations, since the exponent varies at the scale 1/H, which characterizes the off-diagonal elements. 
One may wonder what happened to these intrinsically quantum oscillations in our solution (j4.14|) 
and ([4.17)1 . In order to answer this question we consider the solution to (jAljl 



iA u (k : x) 



d x'G r e (k;x,x')Si2(k,x') 



k 



dtfSxi (k, x - — (t - t'), if — exp - 



i S(Mj 



2 k 



(t-0). (A4) 



The shift in the position, 5x = v(t — t'), corresponds to the retardation for a particle moving with the 
velocity v = k/ko- A similar retardation shift results when the ^-derivative is included in (|A1IIA2)) . 
Since <Si2 is varying very slowly when compared with the oscillatory term, we can expand around 
if = t, 

Si 2 (k,x-^{t- t'\if) = S 12 {k, x) + {if -t)(d t + ^- V)s 12 (k, x) + 0(d 2 ) 
and integrate ()A4j) to obtain 

iA 12 (k,x) = 

This converges provided 



2i 



-S 12 (k,x) 



5{Ml 



2\2 



k-dS 12 (k,x) + 0(d 2 ). 



(A5) 



(A6) 



k ■ d < 5{M ( 



(A7) 



which is precisely the criterion for the validity of the gradient expansion in the off-diagonal equations. 

One may argue that the derivatives acting on the mass must be included as well. To study the role 
of these derivatives, we assume that both the source and the masses depend on the time variable 
only. In this case the off-diagonal equation of motion reads 



[hdt + \{d t M 2 d )d k0 + l -6(M 2 ) - ik 5(E t )]iA 12 (k , t) = S 12 (k ,t) 



(A8) 
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for which the retarded Green function equation is 

[k d t + ~(9tM 2 d )d k0 + l -5{M 2 d ) - ik 6(E t ) + e\G r e (k , k' ;t,t') = 5(k - k' )S(t - t') . (A9) 
When written in the Fourier space 

G:(k ,k' ;t,t') = 

equation (JA9j) becomes 



(Ik 
2^ 



(A10) 



[ ikoK + -(d t M 2 )d ko + U(M 2 d ) - ik 8(E t ) + e]g r e (ko, k' ; k) = 5(k - k' Q ) 



(All) 



% 

2 v-*~- 0/--XJ ■ 2^ 

where we have ignored the higher order gradients in time (that is we took d t M d and 5{M d ) to be time 
independent). The retarded solution of ()A11|) is given by (the solution proportional to —9(k' Q — k ) 
corresponds to the advanced Green function), 



g r e (k ,k' ]K) 



d t M; 



: 6(k - k' ) exp 



d t M 2 



(K-S(E t ))(k' 2 -k 2 ) + (S(M 2 )-2ie)(k' -k ) . (A12) 



This solution is not unique. Indeed, replacing k' 2 by k' [ak' + (1 — a)k } is also a solution, which 
however differs from (jA12|) at higher order in gradients, and hence this ambiguity is irrelevant. 
From ()A12)) we then have 



9tM d 



G r e (k ,k' ;t,t') 



such that Eq. (|A8|) is solved by 



9 / p _ u * 

Z —6{k -k' Q )5(t-t'-^^ 



dtM 2 



> [ J° ] \ 5{M 2 ) - % + k )5(E t ) - 2, e ] 



(A13) 



iAu(k, x) 



dk' dt'G r e (k , k' ; t, t')Si2(k' , i') 
k ° ' b 2 — h' 2 \ 9 



(A14) 



dk' S 12 [k' 0} t 



d t Mj Jd t M d 



d t M3 



The shift in time 8t = —(k 2 — k' 2 )/(d t M}) corresponds to the time retardation caused by the force 
term. Since in the limit of a very slowly varying field (d t — > 1/T — > 0), S12 is varying very slowly in 
comparison with the rapidly oscillating term, we can expand around if = t and k' = ko, 



^ 2 { k '^- dtM j 

and integrate (jA14|) to obtain 

2i 



2 u> 2 - / k' 4- kn 

= S 12 (k ,t) + (k' -k )(d k0 + " 



k 2 — k 



dtM 2 



d t )S 12 (k 0l t) + O(d 2 ,d 2 



t ' u ko 1 



(A15) 



iA 12 (A;o,t) 



5(Mj) 



1 + 



%dt + \(d t M 2 d )d k0 -ifco5(St))]5 12 (fc ,t) + O(0 2 ). (A16) 



One of the higher order terms we have dropped is, for example, 4i[(d t Mj)/(SMj) 3 }d t S(k ,t). This 
converges provided 

k d t , {d t M 2 )d ko « 8(M 2 ). (A17) 

From this and Eq. (jA6HA7|) we conclude that the criterion for validity of the gradient expansion in 
the off-diagonal equations reads 

k ■ d, (dM 2 ) -d k ,k- 5(E) < 5{M 2 d ). (A18) 
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